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On the diffusion of radiation from a point or a line source 
in an infinite medium 
J. T. JErFrERIES 
Division of Physics, National Standards Laboratory, Sydney. 
SuMMARY. — The approximate equation of transfer for diffused radiation in a scattering medium is oblained and applied to the 


cases of infinite media illuminated by a point source or by an infinite line source of radiation. 


SOMMAIRE. — Une relation approchée est obtenue pour la propagation du rayonnement diffus dans un milieu diffusant, puis appli- 
quee aux cas de milieux infiniment étendus éclairés par une source ponctuelle ou par une source linéaire de longueur illimitée. 


ZUSAMMENFASSUNG. — Die Naherungsgleichung fiir die Ausbreitung der diffusen Strahlung in einem streuenden Medium wird 
aufgestellt und auf den Fall eines unendlich ausgedehnten Mediums angewendet, wobei die Strahlungsquelle entweder als punkt- 


formig oder als linear unendlich angenommen ist. 


I. The equation of transfer. — CHANDRASEKHAR [1] 
has shown that the intensity of radiation in an iso- 
tropically scattering medium is given by the equation 
of transfer 


where J = I(x, y, z; 1, m, n) is the specific intensity, 
x the extinction coefficient, mw) the albedo for single 
scattering, and /, m, n are direction cosines. The 
source function P is defined as the ratio of the rate of 
emission of radiation per unit volume per unit solid 
angle to the extinction coefficient x. In a purely scat- 
tering medium, the source function P arises from 
scattering of direct radiation from sources inside and 
outside the medium. 

The exact solution of (1) would be very difficult ; 
however, by using a method similar to that used by 
Eppincron for plane parallel media, a good approxi. 
mation may be obtained for three dimensional pro- 
blems as follows. 

Integrating (1) over a solid angle 4 = it follows that 
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where Fx ==) 1) dw is the « component of the flux 
of radiation, and F and F are defined similarly. 
nS fey dw (all integrals are taken over the whole 
sphere). Now substituting for / from (1) it 1s found 
that, since [1 do = 0, 
oe = A) a) 
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which may be written 
1 [oKar | OKay | Kaz 
x ee | a] Fe ; 


The approximation introduced by Epprneron and 
adopted here consists in taking 


Ku = Ky= Ka = 13 [ Ido = Js; 
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here, however, there is the complication of ‘“‘ mixed ” 
integrals of the form Kx, which, although small for J 
roughly isotropic, may be comparable with those of 
the form Kzz in many cases of interest such e. g. as 
near a point source. If we restrict attention to the 
diffused component, however, the. radiation field 
should be reasonably isotropic ; and so, to a good 
approximation, Kay = Ky: = Ka = 0. 
Substituting in (2) it follows that, if A = 1 — mp, 
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which is the analogue of the MILNE-EDDINGTON equa- 
tion for a three dimensional medium. 

As examples of the application of (4), the cases of 
point source and line source illumination in an infinite 
medium are considered. 


II. Point source illumination. — Consider a point 
source emitting energy at a rate 4x Q. At a distance r 
from the source, 
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so that (4) becomes, in terms of the coordinate r, 
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To solve (6) we apply the method of variation of 
parameters (ForsyrH [2]) and try a solution of the 
form 


(7) J(r) = A(r) 


where k = oe and A(r) and B(r) are determined’ 
in the usual way from the equations : 
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where C = 3x? mQ. 
Solving for A and B, the solution (7) becomes 


@ e Ar or ee (%—A)r 
"iC fe | dr 
(9) po eee 
a el” Wake e (*+h)r 
a ee ——— dr. 
2, k ES ay Ife 


The integration constants a, and a, are to be deter- 
mined from the boundary conditions that, as r >a, 
J > 0 and that as r—> 0 the radial flux of the diffused 
radiation across a spherical surface surrounding the 
source is given by the equation 
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The convergence at infinity clearly requires that 
d, = «©. The second condition gives, making use of 
the series representations for the exponential inte- 
gral : 
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(10) Ei [— (xk) a,].= log ane 
x—hk) 
where 
(11) HA(— 2) =i du. 


The solution (9) appropriate to this problem then 
takes the form 
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Numerical values of the solution (42) are shown in 
Table 4 for 2 = 4/42, 1/75 and 41/300. 

At distances such that xr > 1, equation (12) takes 
the form 
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For high scattering (& small), this becomes 
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J(r) 
which is the solution of equation (4), with P equated 


to zero. This limiting form of the solution for a point 
source in an infinite medium is well known. 


TABLE 1 


Point Source Illumination 


ley KAA, De — A) A = 1/300 
0 ea) a) ora) 
0.1 1.728 1.168 0.763 0 
0.2 LSS 0.678 0 0.391 5 
0.4 0.614 5 0.3188 0.167 1 
0.6 0.392 8 0.1818 0.091 6 
0.8 0.260 8 0.113 0 0.056 2 
1.0 0.184 5 0.074 4 0.036 9 
? 0.122 0 0.050 8 0.025 2 
1.4 0.090 8 0.035 7 0.017 7 
1:6 0.066 1 0.025 6 0.012 7 
1.8 0.048 7 0.018 6 0.009 2 
2.0 0.036 2 0.013 8 0.006 8 


Values of J/3 x? Q a as a function of kr for various values 
of 2. 


JII. Line source illumination. — Consider a 
‘* transparent ”’ line source emitting energy at a rate 
4xQ per unit length and take the source as the z axis 
of a cylindrical coordinate system. The quantity P 
will be a function of e only and may be calculated as 
follows. 
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The radiation emitted per unit solid angle per second 
from unit volume at the point S (fig. 1) due to scatter- 
ing of radiation from dz is 
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where c = w, x is the scattering coefficient. The func- 
tion P(g) is then given by 
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or, in an equivalent form, 
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The transfer equation (4) then becomes, in cylin- 
drical coordinates : 
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withk = V3 axandC'’ = 620m. 
The solution of this equation may be found in a 
similar manner to that used above ; thus 
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K, and J, being the BEsseL functions of purely ima- 
ginary argument and order zero. 

The boundary conditions are that as ep > », J>0 
and, as e — 0, the radial flux 
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These are easily found to require respectively that 
ad, = om and a, = 0 so that the solution (18) appro- 
priate to this problem becomes : 
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Numerical values of the solution (19) are shown in 
able). fora 4420 1/75, 4j300- 
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TABLE 2 


Line Source Illumination 


ke A= 1/12 A= 1/75 = 1/300 
0 lS 7ASpe Pay! 3.148 
0.1 1.520 2.010 2.306 
0.2 1.304 1.595 1.735 
0.4 0.970 1.065 ile stl 1S 
0.6 0.733 0.755 0.781 
0.8 0.560 0.560 0.567 
1.0 0.430 0.427 0.422 
i 0.335 0.323 0.319 
1.4 0.260 0.247 0.245 
1.6 0.205 0.190 0.189 
1.8 0.160 0.148 0.146 
2.0 0.125 0.116 0.114 


Values of J/6% Qay as a function of ko for various values 
of i. 


For e large, (19) takes the form 


6+ 14 

—1 
which has an interesting resemblance to the corres- 
ponding solution (13) for a point source of radiation. 


For high scattering, @ is large and (2) assumes the well 
known limit. 


(21) J ~ 6% Q a Ko(k 9) 


which is the solution of (1) with P equated to zero 
for a line source in an infinite medium. 


(20) J ~ (6 x Q cy 6/2) log 


| K,(kp) 


Acknowledgement. — The author is indebted to 
Dr R. G. GiovaneLit for helpful criticism and cla- 
rification of certain points in the argument. 


Note added in proof. — It has been drawn to the 
attention of the author that a similar problem to that 
considered here has been discussed by R. W. Ditcu- 
BURN (Proc. Roy. Irish Acad. 44 A, 1933, 47) using a 
SCHUSTER-SCHWARZSCHILD approximation for the so- 
lution of the transfer equation. 
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Transmissions des fréquences spatiales 
dans les systémes récepteurs d’ondes courtes 


Jacques ARSAC, 
Laboratoire de Physique de Ecole Normale supérieure 


SomMAIRE. — Une élude des procédés d’ebservalion ulilisés en radioastronomie, et particuliérement des réseaux d’antennes, conduil a 
des résultals analogues & ceux de Voptique instrumentale classique. Si Von décompose la fonction a étudier en fréquences spatiales, 
on montre qu’une antenne opére comme un fillre pour ces fréquences, et peut étre caraclérisée par une bande passante. La limita- 
lion de celle-ci fait que le résultat d’une observation n'est pas la fonction a éludier elle-méme, mais un polynome trigonomeétrique 
qui en réalise une approximation. 

L’étude des réseaux établit que ceux-ci Lransmettent aulant d’harmoniques distincls qu’il y a de distances distinctes entre an- 
tennes du réseau prises deux a deux de toutes les maniéres possibles. Ceci conduit a chercher des réseaux d’antennes, que nous 
appelons incomplets, formés dun nombre donné d’antennes, et transmettant le maximum possible @harmoniques. Avec 4 antennes 
on transmet 6 harmoniques. Le résultat de Vobservation est alors la somme de FouriER d’ordre 6 de la fonclion observée. De tous 
les réseaux transmettanl 6 harmoniques, c’est celui-la qui fournil la meilleure approximalion quadratique moyenne. Avec un 
nombre plus grand d’antennes, le résultal obtenu n’est plus une somme de FouRiER, mais on peul l’y ramener par une transfor- 
mation linéaire. L’utilisation de tels réseaux permet de lirer le meilleur parti des antennes dont on dispose, car, sans diminuer 
la luminosité, elle permet d’augmenter le nombre d@harmoniques transmis (par exemple 286 au lieu de 32 avec 32 antennes). 
Inversement elle permet de réduire au minimum le nombre d’antennes nécessaires a Vobtention dun nombre donné d’ harmo 
niques (perte de luminosilé mais économie de matériel). 

La considération de certains réseaux incomplets permet de montrer que Von peut améliorer le pouvoir séparateur de cerlaines 
antlennes en les faisant interférer avec dauires plus petites convenablement disposées. 

On montre en conclusion que la notion de lobe secondvire n'est pas la plus importante dans la réalisation @une antenne. 


ZUSAMMENFASSUNG. — Hine Untersuchung der Beobachtungsverfahren, die in der Radioastronomie tiblich sind und insbesondere 
von Antennengiltern fiithrl auf Ergebnisse, die denen der klassischen Instrumentalopltik analog sind, Die FourtER-Zerlegung der 
zu untersuchenden Funktion zeigt, dass die Antenne ftir diese Frequenzen wie ein Filter wirkt, und sie kann durch eine Band- 
breile gekennzeichnel werden. Wegen der endlichen Ausdehnung dieser Bandbreile lieferl das Beobachiungsergebnis nicht die 
gesuchte Funktion selbst sondern ein trigonometrisches Polynom das lediglich eine Annéherung darstellt. 

Eine genauere Untersuchung der Gitter zeigl, dass sie ebensoviele verschiedene FourteR-Komponenten tibertragen wie es 
unterschiedliche Abstdnde zwischen den Antennen des Gitters gibt, wenn man sie auf jede migliche Arl paarweise zusammenfasst. 
Das fithrt dazu, Antennengitter zu suchen — wir nennen sie unvollstdndige Gittter — die aus einer gegebenen Zahl von Antennen 
aufgebaut sind und méglichst viele FourtER-Komponenten liefern. Mit 4 Antennen bekomml man deren sechs. Das Ergebnis der 
Beobachtung ist also die Fourter-Summe sechster Ordnung der Beobachtungsfunktion. Es ist von allen Gillern, die 6 Ober- 
schwingungen liefern, grade dasjenige, welches die beste mittlere quadratische Anndherung ergibt. Mit einer grésseren Zahl von 
Antennen erhdlt man nicht mehr eine FourteER-Summe, ober man kann das Ergebnis durch eine lineare Transformation darauf 
zurtckfithren. Die Anwendung solcher Gitter erlaubl, einen méglichst vorteilhaften Gebrauch ven den vorhundenen Antennen zu 
machen. Denn man kann, ohne die Helligkeit herabzusetzen, die Zahl der itbertragenen harmonischen Komponenten erhéhen 
(z. B. 286 anstelle von 32 mit 32 Antennen). Umgekehrt lasst sich die Zahl der Antennen, die nétig ist, um eine geyebene Anzahl 
von FourteER-Komponenten zu bekommen, auf ein Minimum herabseltzen (Lichtverlust aber Materialersparnis). 

Die Betrachtung einiger unvollstandiger Giller ldsst erkennen, dass die Trennschdrfe gewisser Antennenanordnungen ver bessert 
werden kann, wenn man sie mit anderen, kleineren Anordnungen geeigneler Verteilung zusammenwirken lésst. 


_ Schliesslich zeigt sich, dass bei der Herstellung einer Antenne das sekunddre Maximum nicht von ausschlaggebender Bedeutung. 
isl. 


Summary. — A study of the observation procedures used in radioastronomy, and particularly arrays of aerials, leads to results analo- 
gous lo those of classical instrumental optics. If the function under investigation be decomposed into spacial frequencies, it is shown 
that an aerial acts as a filter for these frequencies, and has a characteristic band-pass. The limitation of this leads to the fact 
that the result of an observation ts not the actual function investigated, but a trigonometrical polynomial approximation to it 

Tt is shown that an array of aerials transmits as many distinet harmonics as there are distinct distances between the aerials 
of the array taken two at a time in all possible ways. This leads to a search for aerial arrays, which we call incomplete, formed of 
a given number of aerials and transmitting the maximum number of harmonics. With 4 aerials 6 harmonies are transmitted 
The result of the observation is then the sixth order FourtEer sum of the observed function. Of all arrays transmitting 6 harmonics, 
this gives the best meansquare approximation. With a greater number of aerials the resull obtained is no longer a Fourier sum, 
bul it can be given this form by a linear transformation. The use of such arrays permits the best use to be made of the aerials 
available, for, without loss of intensity, it allows the number of harmonics transmitted to be increased (e. g. 286 in place of 32 
with 32 aerials). Inversely a given number of harmonics may be obtained with a reduced number of aerials (loss of intensity 
but economy of material). 2 

A consideration of certain incomplete arrays shows that the resolving power of certain aerials can be 
interfere with other smaller ones suitably arranged. 

It is shown that the concept,of secondary lobe is not of primary importance in an earial. 


increased by making them 


Un des problémes importants de la radioastronomie dans un petit angle solide , on a la possibilité de 
est observation de la répartition de brillance radio- _mesurer l’énergie émise par une petite portion de la 
électrique sur la surface de la sphére céleste. Elle se sphére céleste, c’est-a-dire la brillance de cette por- 
fait a l'aide d’antennes qui captent lénergie électro- tion. Or, si on veut que cet angle solide ait pour 
magnétique d’une certaine fréquence et de récepteurs demi-angle au sommet 0, il faut que Pantenne ait des 
qui mesurent l’énergie captée. Si lon s’arrange pour dimensions de l’ordre de a = 1/8 @autant plus grandes 
que Pantenne ne transmette au récepteur, de toute que d est plus grand et 9 plus petit. C’est dans ce but 
Pénergie qui lui parvient, que celle qui est contenue qu’ont été réalisés les miroirs paraboliques géants 
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15 m de diamétre pour 4 = 1 cm en Amérique, 25 m 


pour 4 = 20 cm en Hollande ; 75 m pour 2 = 20 em 
en construction en Angleterre. 
Mais la difliculté de réalisation de ces grandes 


antennes est telle que Pon renonce souvent a obtenir 
une antenne directive dans toutes les directions. On 
se contente dune antenne directive dans un seul plan. 
Une seule de ses directions est tres grande. On mesure 
non l’énergie émise par une petite surface, mais celle 
@une petite bande perpendiculaire au plan de direc- 
tivité de Pantenne [1]. 

La réalisation en étant encore délicate, on fait appel 
aux méthodes interférentielles : 

Méthode de Mrcuetson a 2 antennes (proposée par 
Pawsey et utilisée surtout par Rye [2, 3]). 

Interférometres a plusieurs antennes, tel le réseau 
utilisé par CurisTIANSEN [4]. Ces réseaux ne valent 
que pour létude d’astres brillants sur fond sombre, 
car la présence de plusieurs lobes principaux n’est 
acceptable que si un seul d’entre eux recoit de l’éner- 
gie : le diametre apparent de l’astre doit étre inférieur 
a la distance de 2 lobes. 

De toute maniére il mest pas possible d’obtenir 
que lantenne ne recoive d’énergie que d’une portion 
limitée de la sphere céleste ; elle en recoit de tout 
Pespace, mais la transmet au récepteur avec une inten- 
sité variable. On estime que l’antenne la meilleure est 
celle pour laquelle le lobe principal est le plus étroit 
possible et les lobes secondaires les plus faibles pos- 
sible. DoLpH a calculé un type d’antennes optimum [5] 
a largeur de lobe principal donné, il a les lobes secon- 
daires les plus petits possible, a niveau de lobe secon- 
daire donné, il a Je lobe principal le plus étroit possible. 
Nous ne pensons pas toutefois que cette maniére de 
concevoir le probleme fournisse les meilleures solutions. 


Equations du probleme. — Rappelons la définition 
du gain d’une antenne : une onde électromagnétique 
plane de flux d’énergie unité arrivant sur l’antenne 
dans une direction repérable par 2 variables 0, 9, dans 
des axes liés a l’antenne, celle-ci transmet au récepteur 
une énergie G(6, 9). Si ’antenne n’est directive que 
dans un seul plan, G(9, ¢) est indépendant de 9. 
Nous appellerons « gain de Vantenne » la fonction 
Bo Gi 8), 

_ Avec une telle antenne, on ne peut mesurer la bril- 
lance en chaque point, on ne peut connaitre que la 


fonction 
f(0) = [ B(0,9) de. 


Nous prendrons pour variable 9 l’écart angulaire 
sur un grand cercle de la sphere céleste 6 étant la 
position de l’axe de l’antenne, l’énergie /(~) dy émise 
par une tranche de precue dy a labscisse v est regue 
avec un gain G(— 6 + 4) = g(0 — ¥). Au total : 


—i. 
Crest l’équation de la formation d’une image par 


TRANSMISSIONS DES FREQUENCES SPATIALES 113 


un instrument d’optique. L’interprétation d’une telle 
équation est classique [6]. Si Von appelle 7(f) la 
transformée de Fourier de la fonction f, la relation (a) 
equivaut a 


(b) IE) (2) 


De plus si A(x) est la répartition d’amplitude pupil- 
laire 
fen atin | A(x). A(w — u) da. 


Nous ne pouvons utiliser sans précautions de tels 
résultats. Dans le cas qui nous intéresse, les fonctions 


f et g sont périodiques de période 2 x. La transformée 
de Fours de Pune d’elle 
(+o re 
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est nulle si w= p/2 x et infinie si w = p/2 x. On ne peut 
la considérer comme une véritable fonction. I] faudrait 
faire intervenir la théorie des distributions [7]. 

Etant donnée la périodicité de f et g, il y a leu de 
les considérer comme des séries de FourreEer. 
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%, apparait comme la valeur pour Vargument 
u == p/27 de la transformée de Fourrer de la fonc- 
tion égale af si| 8| <x, nulle si| 6| >7. 

Il est done inutile de faire intervenir une fonction 
T(f), une série discrete de valeurs suffit. Il en est 
de méme pour g(9). 

Mais nous ne pourrons pas en général considérer les 
ad, comme les valeurs pour certains arguments de la 
fonction définie par (c). La théorie des antennes donne 
en effet : 

2 riz sin 0 ° 


I(z)e 4% dx 
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Sauf pour des antennes telles que l’on puisse négliger 
g(0) pour des valeurs méme petites de 6 (ce sera 
moyennant certaines précautions le cas des réseaux) 


nous ne pouvons assimiler 9 et sin @ dans cette équa- 


WAS 0 
tion sans introduire d’erreur. Si l’on développe e. 7" 


en série de FOURIER 


+2 
iZsin 0 a , ik 0 
Ce eG CAR aaa 
il vient 
fa a *(y Of = ap 
Je Caan =) dx dy. 


Une ee importante de ceci est que les 
coefficients a, deviennent négligeables a partir dun 
certain rang, mais non rigoureusement nuls [comme 
on Vaurait déduit de (c¢ I. Crest pour cela qu il est 
possible de réaliser des antennes superdirectives pour 
lesquelles la largeur du lobe principal est inférieure 
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ala valeur A/a. On sait [8] que de telles antennes ont 
un tel facteur de surtension qu’elles sont inutilisables 
en pratique. 

De toute maniére, la précision des mesures étant 
limitée, on ne peut connaitre les valeurs de a, au dela 
d’un certain rang, car elles deviennent de l’ordre des 
erreurs de mesure. La suite des a, est done limitée. 
Au total : 
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1 ip 0 
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Une antenne apparait ainsi comme un filtre de 
fréquences d’espace, transmettant un certain nombre 
dharmoniques de la fonction /(0), en modifiant leur 
amplitude et leur phase. La suite des a, Jone pour une 
antenne le méme role que la bande passante pour 
un amplificateur. Nous Vappellerons « bande pas- 
sante spatiale » de Pantenne, 

Considérons un astre de diametre apparent inférieur 
ou au plus égal a 2 6 sur fond sombre. La fonetion 
/(®) correspondante est nulle si} 0|> 6. On pourra la 
representer par 


= ©9 ipt 0 
ms se laa 
(8) = Yi ane 6 gi |s0:} 06 
— 6 i> 
= 0 s1_ | 8 | > 6 
Les seuls harmoniques nécessaires a la connaissance 


Py, ; 
de /(®) sont ceux de rang — .Or unréseau d’antennes 


y) 
omnidirectionnelles, de gain Ep a des abscisses multi- 
ples entiéres dune méme longueur fondamentale a 
%, = E£, @a pour gain 


Kk i | Ky a sin 8 |2 
~ 

BG) => noes na 
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Pour 6 petit, on peut confondre g,(0) avec 


K 2ri LE, a8 2 

v = ae 
g(6) = 24) Gs e X 
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g(8) est périodique de période A/a. Si a/a est assez 
grand, on pourra confondre @ et sin 9 pour| 6|<4/2a, 
c’est-a-dire sur une période de g(9). On peut donc 
substituer a la fonction gp(9) la fonction g(@) définie 
comme suit : 

10] @|<r/2a, g = gr. 

2° g(0) est périodique de période 2/a. 

g est un polynéme trigonométrique de rang WV 


+N 2Ti pad 
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Si a/A = 1/2 B les harmoniques transmis par le réseau 
sont ceux suffisant 4 la détermination de /(@). En ré- 
sume : 
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Etude de l’approximation réalisée. Nous ne pouvons 
connaitre f(0). Nous Vapprochons par P(®), poly- 
nome trigonométrique d’ordre N. P(6) réalise un 
procédé de sommation de /(@) [9]. On sait que le 
procédé de Fourter fournit la meilleure approxima- 
tion quadratique moyenne mais non forcément la 
meilleure approximation en tout point. TOEPLITZ a 
introduit de maniére systématique les procédés de 
sommation : si Pon a 


—- NN 
Fe EN NN ik@ 
n(y, 9) VS USIN a © 
N 
sous des conditions tres générales lim, 
: IN > co 


Tn(y, 9) = f(®) 
Les résultats de ’analyse peuvent se résumer ainsi: 
e étant le maximum de | 7 — f |sur Vintervalle 2 x, 
on peut trouver un procédé de sommation rendant 
¢ minimum. Mais ce procédé dépend de la fonction 
(6). D’autre part le procédé de Fourter est le seul 
reproduisant sans erreur un polyndéme trigonométrique 
@ordre inférieur a N + 1. Ces résultats sont rappelés 
et mis en valeur dans un récent travail de J. C. 
Simon [10]. Devant Pimpossibilité pour nous d’adop- 
ter un critere d’erreur locale, nous adopterons le cri- 
tere de Perreur quadratique moyenne. 
On rendra minimum cette erreur par le choix de la 
somme de Fourrer. a, = 1. La fonction g(@) résul- 
te aN 
sin ——__—. 9 


tante est g(@) = z 


Comme elle peut 
sin > 


prendre des valeurs négatives, elle est inacceptable. 
On leve la difficulté en acceptant que P et f puissent 
différer d’une constante additive et d’une constante 
multiplicative (on se contente de mesures relatives 
ou lon fait un étalonnage). On peut alors calculer 
Pécart quadratique moyen entre f et A + BP, et 
choisir A et B de maniére a le rendre minimum. Le 
résultat est le suivant : 


; C Lp | ot, 


ay p20 


ie een See 

| : ¥ (x, +H) | %F IPI<N 
E’ = |, ? ilp|>N 
avec @,=%)+1Y, 
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E’ est un écart inévitable di au nombre limité d’har- 
moniques transmis. On constate que l’on diminue 
toujours Perreur supplémentaire due a la forme de 
la bande passante spatiale en prenant Y¥p=O0 cest-a- 
dire g(®) paire. En limitant notre discussion a ce cas 
la, il reste : 
2 (2 ap/ ap /)? 
i s 


a lg 
= or [~,| 


E pea 


nw 


= 2 | ts : 
IPI<N 

Cette expression est nulle si tous les dp, sont égaux. 
Mais a ne figurant pas dans cette expression, on peut 
le choisir de maniére a rendre g(6) constamment 
positive. On peut calculer une antenne répondant a 
ces conditions. Le probleme admet du reste plusieurs 
solutions. 

De toute manieére, il est possible d’éliminer par un 
calcul l'intervention de la forme de la bande passante 
spatiale. En jouant sur le fait que P(6) est un poly- 
nome trigonometrique d’ordre N, on peut calculer a 
partir de la donnée de P(6) les coefficients Bp Gp. SI 
a, est différent de zéro, on peut diviser par a, et recons- 
tituer le polyndme de meilleure approximation qua- 
dratique moyenne. 


; I See Neth 2x7ad 
Si Pon pose / == : A ay Sk eee 
(a) po e (9) 9 NHi = dp a) = ; 
h(g) est un polynéme tel que 
. 2N+14 
FB peed sin 5) o) 
J AY) le — 9) do = - 
cee (2. N + 1) sin > 


Alors le polynéme de meilleure approximation qua- 
dratique moyenne est donné par 


, 2yn if 2q% 
—— WY 9D ———— A 


Le calcul de A(g) se fait une fois pour toutes puisque 
h(e) ne dépend que de la constitution de l’antenne. 
On peut aussi imaginer un dispositif optique faisant 
la méme opération. 


Etude des réseaux. — Nous avons vu que moyen- 
nant certaines précautions la fonction gain du réseau 
N QriyaO |2 
est donnée par g(6) =| © Gye 
1 | 
N 2in pad 
Poposant 70)— >, ee 92% 5 g(0) = 118). (9) = 
1 
27ri (p—q)a9 
sy a 
= 2(0)— » Dd & %& rN 
p 4 


en nous limitant au cas ow les coefficients cp sont réels. 

Par conséquent, deux antennes quelconques du 
réseau fournissent un harmonique de rang propor- 
tionnel a leur distance avec une amplitude égale au 
produit de leur gain. 
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Un réseau transmet autant d’harmoniques distincts 
quil y a de distances différentes entre antennes du 
réseau prises deux a deux, de toutes les maniéres 
possibles. Considérons par exemple le réseau uniforme 
formé d’antennes identiques aux abscisses 0, 1, 2,..., N. 
L’harmonique de rang 1 est réalisé de N —1 manieére- 
différentes (puisque la distance 1 est obtenue entre 
deux antennes consécutives), chaque fois avec l’am- 
plitude 1 (puisque toutes les antennes ont méme gain). 
De méme, ’harmonique de rang pa amplitude N—p. 
Il en résulte que le réseau uniforme forme la somme 
de Féser dordre N de (6). Pour des fonctions trés 
rapidement variables (dans l’espace) le phénomene de 
Gripes est moindre que celui di a la somme de Fov- 
RIER, mais les variations rapides de f(@) sont moins 
bien reproduites. 

On peut se proposer de réaliser un réseau, formé 
d’antennes toutes identiques (dans un but de simpli- 
cité de réalisation : il est tres difficile d’ajuster des 
gains d’antennes dans un rapport donné et cela entraine 
une perte d’énergie) et tel que toutes les distances 
entre antennes soient distinctes. Chaque harmonique 
transmis par le réseau l’est alors une seule fois, et 
comme toutes les antennes ont méme gain, tous les 
harmoniques sont transmis avec la méme amplitude. 
De plus on utilise au mieux les antennes du réseau 
car aucune autre combinaison ne peut permettre d’ob- 
tenir plus d’harmoniques. 

-On est ramené a un probleme d’arithmétique. For- 
mer une suite d’entiers E,, Ey,..., E,, telle que les diffé- 
rences positives de ces entiers deux a deux dj; =E; — E; 
donnent, sans omission ni répétition, la suite des en- 
tiers de 1 aN. 

Avec 4 entiers on obtient la solution 0, 1, 4, 6. Si 
le nombre K d’antennes est supérieur a 4, le probleme 
n’admet pas de solution. Supposons en effet que lon 
ait pu en trouver une. Avec K antennes on peut 


K(K — ae 
obtenir Aiba distances distinctes : 
K(K — 1 
a =1silp|< S ) 
, <(K — 1 
a =Osi|p|> — 
La distance 0 est réalisée K fois : ay = K. 
. K?—K+1 2raeé 
sin 
Alors g(8) = K— 1+ Z E 
Ors (Ge aap. 
sin. ——— 
‘i oN 
pour K > 4 et E “ ies ao cette fonc- 


2 (K?—K +1) 


tion est négative, done le probleme n’admet effective- 
ment pas de solution. 

On peut néanmoins chercher a réaliser des réseaux 
incomplets (dérivant du réseau uniforme par la sup- 
pression d’un certain nombre de ses antennes) afin de 
nutiliser que le nombre d’antennes strictement néces- 
saire a l’obtention d’un nombre donné d’harmoniques. 


3. 


{16 Js 


La bande passante spatiale ne pourra pas étre cons- 
tante. Puisque les différences dij ne peuvent repro- 
K(K — 1 , 

Bees la suite 


duire la suite des entiers de 1 a 


de ces différences présentera des omissions ou des 
répétitions. On ne peut accepter d’omission ; nous 
avons vu que l’on ne pouvait reconstituer le poly- 
néme de meilleure approximation que si les a,, étaient 
différents de 0. 

Le probléme est alors de former une suite d’entiers 
E,, E,, ... Ex telle que les différences positives de ces 
entiers deux a deux reproduisent sans omission, mais 
peut-étre des répétitions, la suite des entiers de 
1 aN, N étant le plus grand possible. Un tel probleme 
n’a pas recu de solution. Nous en avons recherché les 
premieres solutions. On trouve ainsi : 


ue S14 tees N= 6.... 0 répétition 

iO ay aie Ni = 9). Werepetition 

Ine aegis N = 13.... 2-répétitions 
ee ern ee N-= 17.... 4 répeétitions 
| es aie areca N = 23.... 5 répetitions 
Ihe ea aes eee N = 29.... 7 répétitions 
[geet KE een ue ane N = 36.... 9 répétitions. 


On peut obtenir une borne inférieure de N. 

Lacsinte: O14. one Thy 2any 23s eae, 
n+p (n+ 2),....n + (n— 1) (n+ 2), n?4+3n—1 
fournit la suite des différences de 1 An? + 3n —1. 
K? ++4K—9 

I . 

Mais pour K = 10 cette borne donne N = 33 alors 
que nous avons pu obtenir N = 36 ; de méme pour 
Ke=11, elle donne N = 39 alors que 0,4, 3,6; 13, 20, 
27, 34, 38, 32, 43 est une solution possible (nous ne 
pouvons pas affirmer que cette suite fournisse le 
nombre maximum d’harmoniques possible). 

De toute maniere il est certain que le nombre d’har- 
moniques que l’on peut transmettre avec un nombre 
K d’antennes suffisamment grand est de Vordre de 
K2 
Tee Ainsi avec 32 antennes, la répartition uniforme 
utilisée par CHRISTIANSEN lui a permis d’obtenir 
32 harmoniques, alors qu’une disposition différente 
aurait permis d’en obtenir 286 au moins. 

Les réseaux incomplets que nous venons de décrire 
ne fournissent pas la meilleure approximation quadra- 
tique moyenne. Pour en tirer pleinement parti, il faut 
faire subir une transformation linéaire a P(0). 

Il faut considérer la méthode comme un compromis 
entre Vinterférométrie a 2 antennes et l’observation 
par les grandes antennes. 

Comme Jlinterférométrie, elle fournit un certain 
nombre d’harmoniques de la fonction {(0), mais elle 
les fournit tous dans un temps tres court (le temps 
de passage de l’astre devant un des lobes du réseau). 
Elle permet étude de fonctions relativement rapi- 
dement variables dans le temps, car il suffit qu’elle 
soit constante pendant le temps de passage. La mé- 


On peut done affirmer que N > 
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thode est utilisable pour le Soleil en présence de taches, 
mais non avec des « bursts ». 

La réalisation de la mise en phase des antennes d’un 
réseau incomplet étant la méme que pour un réseau 
complet, on peut réaliser une telle opération. On est 
str alors des phases des différents harmoniques de la 
fonction /(0), chose qu’il est pratiquement impossible 
dobtenir avec Vinterférométrie. 

Le résultat P(6) est une approximation souvent 
acceptable de f(0). Si ’on veut corriger Vintervention 
de la bande passante spatiale on est conduit a des 
calculs qui n’exigent pas plus d’opérations que le 
dépouillement de linterférométrie. Ces calculs sont 
utiles méme avec une antenne continue (miroirs ou 
réseaux uniformes) si l’on veut en tirer le meilleur 
parti. 

Si l’on se donne a priori un certain nombre d’an- 
tennes, le réseau irrégulier transmet un plus grand 
nombre d’harmoniques que le réseau uniforme. Une 
disposition irréguliére permet done Vobtention, avec 
le méme matériel et dans les mémes conditions (par 
exemple méme distance entre lobes principaux) d’un 
meilleur pouvoir séparateur. La « luminosité » est la 
méme que pour le réseau uniforme puisqu’elle ne 
dépend que du nombre d’antennes. Si Pon dispose 
a priori d’un nombre donné d’antennes et d’une place 
donnée, la disposition irréguliere permet de choisir 
pour le réseau une longueur fondamentale a plus 
petite. Elle permet Vexploration de zones plus com- 
plexes et plus étendues de la sphere céleste. 


Le réseau 0, 1, 4, 6. — I] convient de faire une 
place a part a ce réseau. I] transmet en effet avec une 
amplitude égale les 6 premiers harmoniques de /(6), 
Crest le seul réseau incomplet parfait. Du point de 
vue de lapproximation quadratique moyenne, c’est 
la meilleure des antennes de méme longueur que l’on 
puisse réaliser. Toute transformation sur P(6) en vue 
daméliorer la précision s’avére inutile. 

Nous avons expérimenté un tel réseau pour l’étude 
du Soleil en ondes centimétriques de longueur 
A=3,2 cm. Les antennes du réseau sont des miroirs pa- 
raboliques de 1,10 m de diamétre. Elles sont telles que 
leur directivité est 2° dans le plan méridien et 7° dans 
un plan perpendiculaire. La distance fondamentale 
du réseau est a = 58 2 (185 em). La longueur du réseau 
est 11,10 m. Son pouvoir séparateur est sensiblement 
de 6’. Il est supérieur d’environ 20 % Aa celui d’un 
miroir cylindro-parabolique de méme longueur. 

Nous donnons (fig. 1) la courbe représentative de 
g(®) pour un tel réseau et sur le méme graphique 
celle relative au réseau uniforme a 7 antennes trans- 
mettant le méme nombre d’harmoniques. On cons- 
tate pour le réseau incomplet de trés forts lobes 
secondaires (30 % du lobe principal au lieu de 4,5 % 
pour le réseau uniforme). Le réseau incomplet appa- 
rait ainsi comme une antenne désapodisée par rap- 
port aux antennes classiques. Ceci va a l’encontre des 
recherches classiques a ce sujet qui tendent a une 
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Fie. 1. 


3 Oo = - 

gO) * A 

—réseau 0, 1, 4, 6. 
nreseal Gad 253, 4.50, 


réduction des lobes secondaires méme au prix d’un 
petit élargissement du lobe principal. 

Nous donnons aussi (fig. 2) Venregistrement du 
passage du Soleil devant une telle antenne. Sur la 
partie centrale, nous avons superposé plusieurs pas- 
sages a des jours différents. Vues les tres faibles 
énergies recues, l’intervention du bruit de fond est 
considérable (le signal recu est lui-méme de la nature 
du bruit). 


Solutions mixtes. — Entre le réseau ot lon a réduit 
au minimum le nombre d’antennes nécessaires a l’ob- 
tention d’un certain nombre d’harmoniques et le 
réseau uniforme qui possede le nombre maximum 
d’antennes pour la transmission du méme nombre 


Temps 
Fig. 2. 


d’harmoniques, il y a place pour un grand nombre de 
solutions ot. le nombre d’antennes est intermédiaire 
entre ceux des solutions extrémes. Lorsque lon aug- 
mente le nombre d’antennes en maintenant constant 
le nombre d’harmoniques, le niveau des lobes secon- 
daires décroit progressivement pour tendre vers la 
limite 4,5 % du réseau uniforme. I] semble du reste 
que la décroissance soit beaucoup plus rapide au début. 
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La bande passante spatiale se déforme elle aussi 
progressivement pour tendre vers la bande en « trian- 
gle » du réseau uniforme. 

Certaines de ces solutions mixtes peuvent toutefois 
conduire a des solutions nouvelles. Ainsi pour trans- 
mettre 23 harmoniques il faut 8 antennes (par exem- 
ple 0, 1, 4, 10, 16, 18, 21, 23). Si l’on en utilise 9 on 
obtient la solution 0, 1, 2, 3, 4, 5, 11, 17, 23. Or ce 
réseau dérive du réseau uniforme 0, 1, 2, 3, 4, 5 a 
6 antennes, transmettant 5 harmoniques, par l’adjonc- 
tion de 3 nouvelles antennes qui apportent 18 harmo- 
niques nouveaux. On ne changera pas le résultat en 
remplacant le réseau uniforme par une antenne con- 
tinue de méme longueur. On a ainsi la possibilité 
@améliorer le pouvoir séparateur d’une grande an- 
tenne en la faisant interférer avec des antennes plus 
petites convenablement disposées. De la méme maniére 
le réseau 0, 1, 2,..., 10, 11, 23 fournit les 23 premiers 
harmoniques. I] dérive du réseau uniforme 0, 14, 2,..., 
10, 11 par adjonction d’une antenne qui apporte a 
elle seule 12 harmoniques. On peut améliorer le pou- 
voir séparateur d’une grande antenne en la faisant 
interférer avec une antenne plus petite convenable- 
ment disposée. I] semble que lon puisse beaucoup 
attendre de cette technique en radioastronomie. La 
dimension des grandes antennes réalisables est en effet 
déterminée par des considérations de résistance des 
matériaux, et le plus souvent de prix de revient. On 
a la possibilité, une fois en possession d’une grande 
antenne, d’en amélorer les performances par une 
technique qui ne leur fait pas perdre leur luminosité. 
Nous donnons (fig. 3) un graphique comparant les 
propriétés des deux réseaux cités ci-dessus. L’amélio- 
ration du pouvoir séparateur est tres nette. 


Pe: 2ma 0 
Fic. 3. A= Z 


gq (®) 


g (0) x 


= TESCAUO oo eas Uli ay peers 
~reseawiO) 1) 253,450: 


Dans l’étude qui précéde nous n’avons pas envisagé 
Vintervention d’un certain nombre de facteurs qui 
viennent modifier légerement les conclusions précé- 
dentes. Les antennes du réseau ne sont pas en général 
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omnidirectionnelles ainsi que nous avons suppose. 
Elles ont elles-mémes une certaine bande passante 
spatiale plus ou moins étendue. Soit B(w) ceile-ci 
(nous remplacons la succession des valeurs discrétes 
par une courbe continue, ce qui se fait généralement 
sans ambiguité). Les harmoniques transmis par le 
réseau étant ceux de rang 27 pa/a (avec l’amplitude ap) 


gee Qn 
pa 
la bande passante résultante est Sy a,B (u sais ef 


—N 

La bande passante spatiale d’une des antennes du 
réseau se trouve répétée autour de chaque harmo- 
nique transmis par le réseau. On a ainsi la posibilité 
d’avoir une bande passante spatiale transmettant tous 
les harmoniques jusqu’a un certain rang (et non seule- 
ment des harmoniques plus espacés de rang 27 p a/A) 
si la longueur de la bande d’une des antennes est égale 
a Vintervalle entre deux harmoniques transmis par 
le réseau. 

D’autre part, nous avons supposé au début que la 
distance fondamentale a entre antennes du réseau 
était grande par rapport a la longueur d’onde. S’il 
n’en est pas ainsi, l’assimilation de gp(0) et g(9) que 
nous avons faite nest plus légitime. I] faut opérer 
comme nous l’avons fait au début pour les antennes 
continues : on développe e/7sin ® en série de Fourrer. 
Si Pon considére alors deux réseaux pour lesquels les 
positions relatives des antennes sont les mémes, mais 
qui différent par la valeur de a, on peut passer sim- 
plement de la bande passante de l’un a celle de l’autre. 

Si a, est Pamplitude de Pharmonique de rang 
2xpa/a transmis par le réseau pour lequel a est grand 
devant a, et a‘, celle de Vharmonique de rang p transmis 
par un réseau pour lequel a = 4/2, on a 


N 
= a Jp(0) sae > a Jop(kr). 
AL 


dp 


Il en résulte que la bande passante spatiale d’un réseau 
de dipdles distants de 4/2 est relativement plus éten- 
due. Ainsi le réseau 0, 1, 4, 6 transmet 6 harmo- 
niques (de rang 2p a/d soit une longueur 12 a/A) alors 
que le réseau 0, 4/2, 2 4, 3, transmet 20 harmoniques 
(de rang p soit une longueur 20). La conjugaison des 
deux phénomeénes que nous venons de signaler permet 
la réalisation d’antennes 4 lobe unique par association 
en réseau de petits réseaux de dipdles en 4/2. 


Retour sur la notion de « lobe secondaire ». — Nous 
avons souligné dans ce qui précéde que g(0) est un 
polynéme trigonomeétrique d’ordre N. Il en résulte 
que g(8) peut étre déterminé, de maniére unique, 
par la donnée de 2 N + 1 de ses valeurs. Si ’on prend 
mee 
RCS Sy) 


pour cela les valeurs e( on peut exprimer 


simplement g(@) par 


2qn \ 
ce 8 ea 
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» aN 
= 2 ¥ 27a 70 
A(9)= @ ? ie n = 8 ? 
(2 N + 1) sin-~ 
ge IN| ‘ 
fi 2qr 
= \ oO eC — ore 
s(o)= Zt Sq H he ie 7. 
L’équation fondamentale donne alors 
“ie NI y, 9 \ 
ay vi q 6 
AS 2 BaP pean? 
+N ip70 
avec S(O) == >, oy Coe 
—N 


somme de Fourter d’ordre N de /(6). En définitive 
c’est cette somme seulement que peut nous donner 
lexpérience. On constate que ce n’est pas la fonction 
g(8) qui est importante, mais seulement un certain 
nombre de valeurs de cette fonction. Peu importe 
que g(®) prenne des valeurs importantes si ce n’est 
298 
2N+14- 
Ainsi Pantenne g, = 16, g, =3 si| p |< 6 donne 
P(8) = 13 S(6) + constante 
(ce cas est celui du réseau 0, 1, 4, 6). 
Pourtant la fonction g(6) correspondante est 


pas en un point d’abscisse 


an JBiae 
28 
BSG) ae OPE enemy, 
sin 
s 8 
et présente de forts lobes secondaires. Mais ceux-ci 
298 
se produisent entre les points d’abscisse -—~—— et 
P : IN+4 


n’interviennent pas dans le résultat. 

La technique d’étude que nous venons d’esquisser 
doit étre manié avec précaution. Elle peut néanmoins 
rendre de grands services pour l’évaluation de l’erreur 
possible en un point. Elle met bien en évidence les 


‘valeurs utiles de g(@) et rejette au second plan la 


notion de lobe secondaire qui n’est pas finalement la 
grandeur critique dans la réalisation d’une antenne. 
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COHERENCE DE L’ECLAIRAGE ET CONTRASTE DE L’IMAGE 


Influence de la cohérence de l’éclairage sur le contraste de image 
d’un point noir en présence d’un petit défaut de mise au point. 


SERGE SLANSKY 
Institut d’Optique, Paris. 


SOMMAIRE. — Le contraste de V'image @un point noir sur fond blane est étudié en fonction du rapport p = 
riques du systéme d’éclairage et de Vinstrument d’observation. Pour p= 


traste est, 


fon = b/a des ouvertures numé- 
= 1 si Vinstrument est parfaitement stigmatique le con- 


en premiere approximation, indépendant de p, et égal au contraste en éclairage cohérent ; par contre, si V instrument 


presente un léger défaut de mise au point, le contraste dépend de p et atteint son maximum au voisinage de p = 0,87. 


D ofp 7 7 7 “eS Dd 
Pour p> i, le contraste diminue toujours quand p augmente et pour les trés grandes valeurs de p le contraste tend vers sa 
valeur en éclairage incohérent, qui est environ la moitié du contraste obtenu POUrp — ie 


ZUSAMMENFASSUNG. — Der Kontrast im Bilde eines dunklen Punktes in hellem Felde wird als Funktion des Verhdltnisses p = dla 
der numerischen Beleutchtungsapertur zur Apertur des Beobachtungsinstrumentes, untersucht. Wenn man das Instrument als 
punktférmig abbildend vorraussetzt, dann ist der Kontrast fiir p <1 in erster Naherung unabhdngig von p und gleich dem Kon- 
trast fir kohdrente beleuchtung. Wenn aber das Instrument eine kleine Abweichung von der besten Strahlenvereinigung aufweisl 
so hangt der Kontrast von p ab und erreicht seinen Héchstwert nahe bei p = 0,87. Fiir p > 1 verringert sich der Kontrast immer 
wenn p zunimmtl, und flir sehr grosse Werte von p ndhert sich der Kontrast dem Wert fiir inkohdrente Beleuchtung, der etwa 


halb so gross ist wie der fiir p = 1. 


SuMMARY. — The contrast of the image of a black point on a white background is studied as a function of the ratio p = b/aof the 
numerical aperture of the illuminating system to that of the imaging system. For p <1, for a perfectly stigmatic instrument the 
contrast is, to a first approximation, independent of p and equal to the contrast with coherent illumination ; in the presence of a 
small defect of focus the contrast varies with p, and attains a maximum in the region of p = 0.87. For p > 1, the contrast decreases 
and tends toward the value for incoherent illumination, which is about one half thal obtained for p = 1. 


Le caleul du contraste de image d’un point noir 
sur fond blane conduit, comme on sait, a des résultats 
différents selon que l’éclairage est supposé cohérent 
ou incohérent [1]. Nous nous proposons d’étudier ici 
ce qui se passe dans les cas intermédiaires, ot l’éclat- 
rage est partiellement cohérent. 

La théorie de la cohérence partielle, qui a été étu- 
diée par divers auteurs, notamment ZERNIKE [2], 
Hopkins [3] [4], BLanc-Laprerre et Dumonret [5], 
WotrF [6], montre que le facteur de cohérence de la 
lumiere dépend des dimensions angulaires apparentes 
de la source lumineuse, ou de la pupille de sortie du 
systéme d’éclairage si l’on emploie un condenseur. Si 
le systeéme d’éclairage a une pupille circulaire, on peut 
done prendre comme parametre déterminant l'état 
de cohérence partielle de la lumiére ouverture numé- 
rique du systéme (évaluée du cété de lobjet éclairé). 
Comme d’autre part, pour un mode d’éclairage donné, 
le facteur de cohérence entre deux points éclairés 
dépend des positions relatives de ces points, pour 
apprécier dans quelle mesure l’éclairage fourni par 
un certain dispositif est pratiquement cohérent ou 
incohérent il convient de préciser les caractéristiques 
de lobjet éclairé et de ’instrument employé pour son 
observation. 

Nous considérerons ici le probleme suivant : Soit 
un instrument d’ouverture numérique @ qui forme une 
image d’un objet constitué par un point noir (c’est-a- 
dire une petite tache noire de dimensions nettement 
inférieures au pouvoir séparateur de l’instrument) sur 
fond blane (c’est-a-dire uniformément transparent). 
L’objet étant éclairé a l’aide d’un certain systeme 
d’ouverture ) = pa (les ouvertures a et b sont éva- 
luées toutes deux du cété de l’objet éclairé, et corres- 
pondent done respectivement a la pupille d’entrée de 
Vinstrument d’observation et a la pupille de sortie du 


systeme d’éclairage) nous cherchons comment varie 
le contraste de Pimage obtenue en fonction du rap- 
port p = D/a. 

Nous nous limiterons ici au cas de systemes centrés 
a pupilles circulaires, et nous considérerons plus par- 
ticulierement le cas oti Vinstrument d’observation 
est un instrument parfait et celui ou il présente un 
léger défaut de mise au point. On sait d’autre part 
que les aberrations du condenseur n’ont aucune in- 
fluence [2] [3]. 

Les caleuls seront faits en utilisant approximation 
de la théorie des ondes scalaires. 


Expression de la répartition d’intensité dans l’image 
d’un objet noir sur fond blanc. — Considérons un 
objet plan (©) de transparence complexe E(M) dont 
on forme une image (J) a Vaide d’un instrument 
optique (fig. 1), que nous supposerons étre un sys- 
téme centré dont laxe optique est perpendiculaire au 
plan de Vobjet. Soit F(M’‘) la fonction donnant la 
répartition des amplitudes complexes dans la tache 
de diffraction produite dans le plan image par un 
point lumineux placé au centre du plan objet. Nous 
admettrons que pour un choix convenable des coor- 
données dans les plans objet et image amplitude 
produite en M’ par un point lumineux placé en M est 
de la forme /(M’‘ —M), c’est-a-dire que les aberra- 
tions varient peu en fonction de la position du point 
objet dans le champ. Nous prendrons ici dans le plan 
objet une unité de longueur égale a la longueur d’onde 
de la lumiére employée, supposée monochromatique ; 
dans le plan image, cette unité sera multipliée par le 
grandissement de Vinstrument. On peut définir aussi 
Vinstrument par le facteur de transmission f(X) de sa 
pupille . La variable X est une grandeur vectorielle 
a 2 dimensions qui définit la direction dun rayon 
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allant du centre du plan objet au point considéré 
de la pupille d’entrée du systeme ; on peut linter- 
préter comme la projection sur un plan perpendicu- 
laire a l’'axe optique du systeme d’un vecteur unitaire 
porté par le rayon (pour un rayon faisant un angle u 
avec l’axe, on a| X |= sin w). Pour les valeurs de X 
intérieures au domaine Q défini par louverture de la 
pupille, ((X) = e®, ou W(X) représente le dépha- 
sage entre ’onde sortant de l’instrument et celle que 
donnerait un instrument parfait (nous supposons 
qu’on a établi une correspondance entre les pupilles 
d’entrée et de sortie de telle sorte que deux points 
situés sur un méme rayon passant par le centre de 
Pobjet aient des coordonnées égales). Pour les valeurs 
de X extérieures a ©, c’est-a-dire correspondant a des 
rayons passant en dehors de la pupille, nous aurons 
f(X) = 0. On sait que F(M’) se déduit de f(X) par une 
transformation de Fourier 


a 


(1) F(M’) = | 1x) e2tiXM! ax. 


ce 


Considérons d’autre part une source plane uniforme 
qui éclaire a travers un certain systeme optique défini 
par le facteur de transmission a(X) de sa pupille ¥. 
Ici, X représente la direction d’un rayon allant du 
centre de l’objet éclairé a un point de la pupille de 
sortie du systeme d’éclairage, cette direction étant 
définie de la maniere indiquée plus haut. Nous sup- 
poserons que l’axe optique du systeme d’éclairage 
est le méme que celui de l’instrument d’observation 
et perpendiculaire au plan de la source. 

Soit A(M) la transformée de Fourter de a(X). 


[ Orr. Acta 


Avee un choix convenable de coordonnées dans le 
plan de la source, amplitude produite au point M de 
Pobjet par un point P de la source sera proportion- 
nelle a A(M—P). En éclairant Vobjet avec une 
source ponctuelle placée en P, nous aurons au point M’ 
de image une amplitude proportionnelle a 


(2) U(M',P) = | A(M—P) £(M) F(M’—M) dM. 

Les différents éléments de la source étant supposés 
incohérents, Pintensité produite par une source éten- 
due s’obtient en intégrant lintensité (carré du module 
de amplitude) produite par un point ; si lintensité 
de la source par unité de surface est prise égale a 
Punité, Pobjet éclairé par Vensemble de la source 
donnera dans le plan image une répartition d’intensité 
de la forme 

(3) I(M') = | | U(M’,P) |2 dP 

Nous supposerons la source suffisamment étendue 
pour que Péclairement résultant (dans la région qui 
nous intéresse) ne soit limité que par louverture du 
systeme d’éclairage, de sorte que l’expression (3) a 
sensiblement la méme valeur que si le domaine d’in- 
tégration par rapport a P était illimité. Dans ces con- 
ditions, il résulte des propriétés de la transformation 
de Fourier qu’en posant 


(4) w(M/,X) =| U(M',P) e2"iXP gp 


Pexpression (3) est équivalente a 
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(5) I(M’) a | u(M’, X) |2 aX 


En combinant les formules (2) et (4) nous obtenons, 
en posant P = M+ Q, 


(6). “a(M-X) = 
7.7 i 

21 + n I] 
[A(—Q) e'1XQ aQ { E(M) F(M’—M) e21X™ aM 
e e 
La premiere intégrale représente l’expression com- 
plexe conjuguée de a(X), comme on peut le voir en 
appliquant les formules d’inversion des intégrales de 
Fourier. En désignant autre intégrale par o(X, M’), 
nous voyons que /(M’) peut se mettre sous la forme 


aX, M) |? aX , 


(3) 9(X,M’) = | E(M) F(M’—M) e?™iMX qx | 
Une expression ¢quivalente a été trouvée par Hop- 
KINS [4] en partant d’une certaine définition du fac- 
teur de cohérence qui fait intervenir la transformée 
de Fourier de | a(X) |? . 

Notons que a(X) mintervenant que par le carré de 
son module, la répartition des intensités dans limage 
ne dépend pas des aberrations du condenseur (repré- 
sentables par une variation de phase sur la pupille 
de sortie) ; ce fait avait déja été signalé par ZErR- 
NIKE [2]. 

En tenant compte que la fonction a(X) est nulle 
a lextérieur du domaine ¥ limité par louverture du 
condenseur et qu’elle est de module unité a lintérieur, 
la formule (7) peut encore s’écrire, en limitant le do- 
maine d’intégration, 

+ 
[| 20%M) |? ax 


a al 


(9) I(M’) = 


Appliquons les formules précédentes au cas ou 
Pobjet est constituée par une petite tache noire (opaque) 
sur un fond complétement transparent. La fonction 
E(M) est nulle a lintérieur d’un certain domaine S et 
égale a l’unité a l’extérieur de S (fig. 2). Désignons 
par E'(M) la fonction qui représente Vobjet com- 
‘plémentaire (trou dans un écran opaque) 


(Aq sur > 


E'(M) = 1 — E(M) = , (.) en dehors de S 


(10) 


Soit. o/(X, M’) Vexpression obtenue en remplagant / 
par E’ dans la définition (8) de ¢(X, M’) 


(11) (X, M’) =f (MOM) = aM 
Nous avons 
(12) @(X,M’) = f(X) e@ "XM — 9'(X,M’) . 


Désignons par g(X) la fonction f/(X) | a(X) |?. Nous 


avons 
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E (M) 


13) | AX) = g(X) sur > 
g(X) ==10 en dehors de & . 
Soit G(M) la transformée de Fourter de g(X). 
En tenant compte des relations (11), (12) et (43), on 
peut mettre /(M’) sous la forme 
(14) 9 I(M’) = J — 2 1,(M’) + 1,(M’) 
avec 
(15) s=f | foo leax = | | gow |? ax 


= « 


(16) 1(M’) = 5 [. LF(M’ — M) G*(M! — My + 
+ G(M' —M) F” (M'—M)] dM 
(7) TIM = I. | o/(X,M’) [2 dX . 


e ~_ 


On peut noter que J représente Pintensité obtenue en 
supprimant la tache noire, et /, Vintensité obtenue 
en remplacant la tache noire sur fond blane par une 
tache blanche sur fond noir. 

Les formules précédentes sont valables pour une 
tache noire de dimensions quelconques, et pour un 
instrument dont les aberrations peuvent étre quel- 
conques (pourvu qu’elles soient représentables par 
le facteur de transmission f(X) indépendant de la 
position du point M dans la région utile du champ de 
Vinstrument). 


Cas ot Vobjet noir est tres petit. — Lorsque la tache 
noire est suffisamment petite, par rapport au pouvoir 
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séparateur de linstrument d’observation, pour étre 
assimilable a un point, on peut admettre en premiere 
approximation que les fonctions /(M'—M) et 
G(M’ —M) ne varient pas lorsque, M’ étant fixé, M 
se déplace sur S. En supposant le « point noir » placé 
au centre du plan objet, nous obtenons 


(eS a S{F(M') G*(M') + G(M’) F*(M!)| 


a 


(18) 


(19) (MX) =X) FM = FM) | PO aM, 


FYE Te] 


S désignant l’aire de la tache noire (en prenant A° pour 
unité de surface) et e'(X) étant la transformée de 
Fourier de E'(M). 

Nous supposerons que Vinstrument d’observation 
(objectif) et le dispositif d’éclairage (condenseur) ont 
des pupilles circulaires, dont les ouvertures nume- 
riques seront désignées respectivement par a et b. 
Le domaine ® sera donc défini par X? < a? et S par 
X? < 6%. La fonction f(X) sera nulle pour X? > @?, 
et pour X? <a? elle aura pour module Punité et pour 
argument une certaine fonction lée aux aberrations 
de Pobjectif ; la fonction g(X) sera égale a f(X) si Pon 
a a la fois X?< a? et X? < b?, et s’annulera partout 
ailleurs. Il en résulte que les fonctions /(M‘) et 
G(M’‘) seront identiques si «< b, mais généralement 
distinctes si a> b. 

Examinons successivement les cas a > b et b > 4a. 
Si a> 6 (le condenseur est moins ouvert que Vobjectif) 
nous obtenons immeédiatement, d’apres l’expression 
(15) de J: = 

(20) B= oitike 


D’autre part, la tache noire ayant des dimensions 
petites par rapport au pouvoir séparateur de lobjectif 
qui est de Pordre de 1/a, il en résulte que la fonction 
e'(X), qui est la transformée de Fourrer de E’(M), 
varie peu sur 2, donc aussi sur XY. Dans les for- 
mules (17) et (19) nous pouvons en premiere approxi- 
mation confondre e’(X) avec e(0) = S, dot nous 
obtenons : 
(21) 


Si b >a (le condenseur est plus ouvert que Vobjectif) 
nous obtenons au lieu de (20) expression. 


(22) J Sng. 


D’autre part, lorsque le rapport p = b/a devient nette- 
ment supérieur a Punité, la formule (21) cesse d’étre 
une bonne approximation, car e’(X) peut arriver a 
varier sensiblement dans le domaine Y qui est plus 
étendu que 2. I] faut alors remplacer |’expression (21) 
par expression plus générale 


ihe; = iS? 


F(M')? . 


(23) 1,(M') = | F(M) [2 [| e'(X) |? dX. 
Si le condenseur est beaucoup plus ouvert que l’ob- 
jectif, b > a, de telle sorte que la tache noire qui est 
petite par rapport au pouvoir séparateur de Vinstru- 
ment d’observation est cependant grande par rapport 
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au pouvoir séparateur d’un instrument dont Pouver- 
ture numérique serait égale a celle du systeme d’éclai- 
rage, on obtient a la limite 


(24) T,(M) = S | F(M’) |. 


Cette expression s’obtient en remarquant que si le 
domaine d’intégration est infini (ce qu’on peut ad- 
mettre si le domaine ¥ est assez grand pour que I’in- 
tégrale de | e/(X) |? prise a l’extérieur de » soit négli- 
geable) on peut remplacer Vintégrale de | e’(X) |? dX 
par celle de | £’(M) |? dM, dont la valeur est évidem- 
ment 8. 


Expression du contraste de !’image. — Nous defi- 
nirons le contraste de ’image du point noir par la 
baisse relative d’intensité que la présence de ce point 
noir au centre du plan objet produit au centre du 
plan image. Une telle définition a un sens a condition 
que les aberrations ne soient pas trop importantes 
pour que /(M’‘) ait un minimum assez net au centre, 
L/expression du contraste sera donc 


(25) y = (t/J) [2 1,(0) — 1,(0)] , 


les grandeurs J, /,, /, ayant les valeurs indiquées plus 
haut. Nous distinguerons deux cas selon les valeurs 
relatives de a et b. 


1° Le condenseur est moins ouvert que lobjectif 
(6 < a). En utilisant les formules (18), (20) et (21) 
nous trouvons 
Se) G*(0)+ G(0) F*(0)] — S?|.F(0) 


(Quand Pouverture } devient tres petite, G(0) se con- 
fond avec x 5? et le contraste tend vers 
(27) y¥, =1—|1—SFO), 


expression équivalente a celle qu’on obtiendrait dans 
le cas de Péclairage parfaitement cohérent. 


(26) y= 


2° Le condenseur est plus ouvert que l’objectif (b > a). 
Les fonctions F et G devenant identiques, /,(M’) se 
réduit a | F(M’) |? S. En utilisant les expressions (22) 
et (23), nous obtenons, comme expression du con- 


traste, 
EE [eX 2 ax| . 


On peut facilement en déduire que, pour une ouver- 
ture donnée de Vobjectif, le contraste devient une 
fonction décroissante de louverture du condenseur 
des que cette ouverture dépasse celle de l’objectif. 

Si le rapport p= b/a nest que légerement supé- 
rieur a Punité, on peut considérer que lexpression (21) 
de L, est encore approximativement valable et écrire 


SOE 


Tae 


(28) 


(29) 


| (0) Re pees ee (yas 


mais cette approximation n’est acceptable que dans 
un domaine assez limité. 
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Si le rapport p = b/a devient tres grand, on peut 
prendre pour /, l’expression (24), et le contraste tend 
vers 

(30) 


2S. 


T™ a? 


= 41 FOP, 
expression €équivalente a celle qu’on obtiendrait en 
supposant l’éclairage complétement incohérent. 

Lorsque les ouvertures de l’instrument et du sys- 
teme d’éclairage sont égales (a = b, p = 1), Vexpres- 
sion du contraste devient : 


(31) Y= 


Comme le « point noir » est supposé avoir des dimen- 
sions petites par rapport a 1/a, nous avons x a2S < 1. 
Done la valeur de y donnée par (31) est A peu pres le 
double de la valeur y; donnée par (30). 

Nous voyons done que lorsque les ouvertures de 
Pobjectif et du condenseur sont égales, le contraste de 
Vimage dun point noir est environ le double du con- 
traste qu’on obtiendrait en éclairage incohérent. I] n’est 
done pas exact d’admettre, comme on le fait parfois, 
que l’on obtient un éclairage pratiquement incohé- 
rent pour p = 1. 

Nous donnons ci-dessous (tableau I) les valeurs 
numériques du rapport entre le contraste y calculé 
pour diverses valeurs de p = b/a d’apres la formule (28) 
et le contraste y; obtenu en éclairage incohérent 
d’apres la formule (30), en prenant comme objet noir 
un petit cercle de rayon < tel que ¢ a= 0,08; on a 
na? S = 0,06. L’expression générale de y lorsque 
Pobjet est un cercle de rayon <« < 1/a est (pour b > a) 


S 
‘ 1g | 2 2 76 
(32) y= —_| FO) PU + Jp (2neb) + J{(2ned)]. 
TABLEAU I 

| p | 1 2 3 4 A 40 0 
| 
| viv; -| 1,94 1,78 1,57 1,88 1,25 1,14 1,06 1,00 

Les résultats sont représentés graphiquement 


(fig. 3); la courbe obtenue en utilisant la formule 
approximative (29) est représentée en pointillé : nous 
-yoyons qu’elle donne une approximation assez bonne 
jusqu’a p = 2 ou 3. 

Les résultats obtenus jusqu’ici ne supposent pas 
@hypothése particuliére sur les aberrations de Vins- 
trument, a part que ces aberrations ne doivent pas 
étre trop fortes, de maniére que l’on obtienne encore 
un contraste appréciable avec un objet dont les di- 
mensions sont assez petites non seulement par rapport 
au pouvoir séparateur de l’instrument considéré mais 
aussi par rapport a celui d’un instrument parfait de 
méme ouverture. Nous allons examiner d’une maniere 
plus détaillée les cas particuliers de instrument par- 
fait et de l’instrument présentant un défaut de mise 
au point. 
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Cas de l|’instrument parfait. — Pour un instrument 


dépourvu d’aberrations et parfaitement mis au point 
le facteur de transmission f(X) est égale a lunité sur 
toute la pupille 


(33) f()) =e POUL ok <raas 
On en déduit 

(34) F(0) = @? , 

(35) GO) aor O* == op? 2-0? 


Pour p< 1 (& <a), nous obtenons 
(36) Veen Tene eae Caer 


Le contraste est done indépendant de louverture 
du condenseur tant qu’elle ne dépasse pas celle de 
Pobjectif. Il en résulte que pour un instrument parfait 
le contraste obtenu pour p = 1 est le méme qu’en éclat- 
rage cohérent. Comme Tobjet est supposé trés petit, 
ra? S <1 et on peut négliger en premiere approxi- 
mation le terme en S?, d’ou 


(37) 


Pour p> 1, le contraste diminue. Nous avons déja 
indiqué plus haut Vallure générale de la décroissance 
du contraste pour les grandes ouvertures du conden- 
seur. Pour les trés grandes valeurs de p nous obtenons 
a la limite la valeur du contraste en éclairage incohé- 
rent. 


(38) 


Les formules (37) et (38) correspondent aux expres- 
sions obtenues précédemment par A. Marecnat [1], 
concernant respectivement I’éclairage cohérent et 
Péclairage incohérent. D’apres ce que nous avons vu 
plus haut, le fait que le contraste en éclairage cohé- 
rent est double du contraste en éclairage incohérent 


es 2S) ae OF 


BG 
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dans le cas d’un instrument parfaitement stigmatique 
résulte du fait que, dans ce cas particulier, le contraste 
pour p = 1 est le méme qu’en éclairage coherent. 


Cas d’un petit défaut de mise au point. — Un petit 
défaut de mise au point produit entre les zones cen- 
trales et marginales de onde un déphasage propor- 
tionnel au carré de louverture, donc le facteur de 
transmission /(X) est de la forme 


(Bo = ela x” pour AZ-—< wu’. 


Le coeflicient « étant proportionnel ala distance entre 
le plan de mise au point utilisé et celui correspondant 
ala mise au point exacte. En posant 


(40) Ce == 0, 


nous obtenons, pour p = b/a <1, 


Pour p> 1, on a les mémes expressions que pour 
p = 1. Nous en déduisons lexpression du contraste, 
pour p< 1, 

ra 


(43) y=2 S— [1 — cos — cos (p*p)+-cos (p?—1)o]. 
(Os is 


Le contraste en éclairage cohérent (p —> 0) est 


sin © 
2 


En éclairage incohérent, on obtient 


(44) a eS ata 


En considérant la dérivée par rapport a p de l’expres- 
sion (43) de y, nous trouvons que pour les valeurs de ¢ 
qui sont inférieures a = la valeur de p correspondant 
au maximum de contraste reste assez voisine de 
p = 9,87. Nous donnons ci-aprés (tableau II) les 
valeurs, calculées pour diverses valeurs de p, du rap- 
port entre le contraste y obtenu avec un défaut de 
mise au point o = x/2 (correspondant a la tolérance 
de RayLercu) et le contraste y) obtenu en éclairage 
cohérent avec Jlinstrument parfaitement mis au 
point. Les valeurs indiquées pour les valeurs de p 
supérieures & Punité ont été obtenues en prenant 
comme objet un petit cercle de rayon 0,08 (1/a), qui 
donne y, = 0,12. 


TABLEAU II 


on 
i) 


io’) 


p | 0 0,20 0,50 0,87 1 


the 0,64 0,66 0,74 0,84 0,81 0,51 0,46 0,41 


EGO. 
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Nous donnons également (tableau IIT) les résultats 
dun calcul analogue pour @ = 3x/4 (défaut de mise 
au point correspondant a une fois et demie la tolé- 
rance de RAYLEIGH). 


TABLEAU III 


P | 0 0,20 0,50 0,87 1 5 10 00 


y/Yo-| 9,30 0,33 0,49 0,65 0,61 0,39 0,35 0,30 


Les résultats sont représentés graphiquement (fig. 4). 
Nous donnons d’autre part (fig. 5) les courbes repré- 
sentant en fonction de p le défaut de mise au_ point 
tolérable 8 pour une baisse de contraste de 20% ou 
de 40°. L’unité employée est telle que 5 = 1 corres- 


TT 


pond ag=—. 


Conclusion. — L’étude de influence des conditions 
d’éclairage sur le contraste de image d’un point noir 
montre que l’éclairage cohérent et léclairage incohé- 
rent peuvent effectivement s’interpréter comme deux 
cas limites de Jléclairage partiellement cohérent, 
limites correspondant respectivement a des valeurs 
infiniment petites ou infiniment grandes du rapport 
des ouvertures numériques du systeme d’éclairage et 
de l’instrument d’observation. Pour une ouverture 
donnée de l’instrument il est en principe toujours pos- 
sible de réaliser un systeme d’éclairage d’ouverture 
assez petite pour que l’éclairage soit pratiquement 
cohérent. Par contre, un systeme d’éclairage ayant 
une ouverture assez grande pour que l’éclairage puisse 
étre considéré comme incohérent n’est effectivement 
réalisable que si l’instrument d’observation est lui- 
méme assez peu ouvert. I] en résulte que, par exemple, 
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léclairage sera assez loin d’étre incohérent pour les 
ouvertures généralement utilisées en microscopie ; 
nous avons vu qu’une ouverture du condenseur égale 
a celle de Pobjectif est encore nettement insuffisante 
pour que léclairage puisse étre considéré comme 
incohérent. Par contre, il arrive souvent que l’éclai- 
‘age soit pratiquement incohérent dans le cas de l’ob- 
servation des objets éloignés : par exemple, pour une 
lunette astronomique qui vise une planéte passant 
devant le Soleil, ’ouverture numérique de Vinstru- 
ment, évaluée du coté objet, est de l’ordre de 10-, 
celle du systeme d’éclairage constitué par le Soleil 
est de lordre de 10-2, dot un rapport de Vordre de 
108, pratiquement infini. 

D’autre part, étude de la variation du contraste 
en présence d’un défaut de mise au point montre que, 
dans le cas ot Pouverture du condenseur ne dépasse 
pas celle de Vobjectif (cas de la microscopie par 
exemple), si Pouverture du systeme d’éclairage n’in- 
flue guere sur la valeur du contraste qu’on peut obte- 
nir avec un instrument parfaitement au point, elle 
influe considérablement sur la baisse de contraste 
produite par un défaut de mise au point. Un choix 
convenable des conditions d’éclairage peut donc in- 
fluencer dune maniere notable la tolérance sur la 
précision de mise au point de l’instrument. 


REFERENCES 


[1] A. MarécuaL, Rev. Opt., 82, 1953, p. 649. 

[2] F. ZERNIKE, Physica, 5, 1938, p. 785. 

[3] H. H. Hopxins, Proc. Roy. Soc. A., 208, 1951, p. 263. 

[4] H. H. Hopxins, Proc. Roy. Soc. A, 217, 1953, p. 408. 
BLANC-LAPIERRE et P. DumonteET, Rev. Opt., 34, 
19555) pes 

[6] E. Wotr, Proc. Roy. Soc. A, 225, 1954, p. 96. 


Manuscrit recu le 23 mat 1930. 


126 ET. 


HAVLICEK 


[Orr. Acra 


On a set of classical methods for the computation of optical systems 


F. I. Haviicek 
J. Sreran Institute, Ljubljana- Yougoslavia 


Summary. — Starting from Gaussian oplics, PetzyaL sum and a first criterion for low distortion for the first layout, there 
are foreseen three steps in the calculation of the system. In the first there is combined a SEIDEL approximation for thin lenses 
introducing ARGENTIERT’S coefficients, with an analysis for a thick system following BEREK’s method. Iteration of this two 
methods with slide rule accuracy give the starting point for the next higher stage of correction of the system of thick lenses. In 
the nexl slep a variation method of SEIDEL coefficients leads with parallel trigonometrical checks to the final stage of correction. 
In the last step is introduced a control of the results obtained with the former methods by using an analysis of the wavefront. 


SOMMAIRE. 
sion, trois stades sont prévus dans le calcul dun systéme. 


En partant de UV optique de Gauss pour élablir en premier lieu la somme de Perzvat ef une condilion de faible distor- 


Dans le premier on combine une approximation de SE1pEL pour lentilles minces introduisant les coefficients d’ ARGENTIERI, 
avec Vanalyse d’un systéme mince selon la méthode de BEREK. La répétiiion de ces deux méthodes (avec la précision de la régle 
a calcul), donne le point de départ pour le 2° stade de correction du systéme de lentilles épaisses. Dans ce 2° stade, une méthode 
de variation des coefficients de SErDEL combinée avec des calculs de contréle trigonométriques conduisent au dernier stade de la 
correction. Dans la derniére phase du calcul on introduit un contréle des résultats obltenus par les méthodes précédentes en analy- 


sant la forme de la surface d’onde. 


ZUSAMMENFASSUNG. 


Nachdem man im Rahmen der Gaussschen Optik einen ersten Ansalz mit dem gewiischten Wert der Pe‘z- 


valsumme und fiir eine geringe Verzeichnung gemacht hat, wird empfohlen, die weitere Berechnung eines optischen Systems in 
drei Schrilten vorzunehmen. Zundchst kombinierl man im Seidelschen Gebiet die von ARGENTIERI eingeftihrten Koeffizienten 
fir diinne Linsen mit der analytischen Methode nach BER®Ex fiir dicke Linsen. Die wiederholte Anwendung dieser beiden 
Methoden mit Rechenschiebergenauigkeit bildel dann den Ausgangspunkt fiir die ndchste Stufe in der Korrektion des Systems aus 
dicken Linsen. Bei diesem Schritt ldéuft der Variation der Seidelschen Koeffizienten die trigonometrische Kontrolle parallel. 
Schliesslich werden die Ergebnisse an Hand einer Analyse der Wellenflache gepriift. 


Methods of designing optical systems, which are 
believed to be of general use, have been developed at 
the J. Sreran Institute. 

Starting from Gaussian optics and the PETzvaL 
condition for a first layout, three steps are used in 
calculating an optical system. First the SEIDEL 
approximation for thin lenses (using ARGENTIERI’S 
coefficients [1]) is combined with a thick lens ana- 
lysis using Berek’s method [2]. Iteration of these 
methods to slide rule accuracy gives the starting point 
for the next stage. Here a variational method of SEI- 
DEL coefficients with parallel trigonometrical calcula- 
tions leads to the final stage of correction. Lastly the 
results obtained by the former methods are checked 
by analysis of the wavefront. 

A first layout should give a low value for the contri- 
bution of the separate surfaces, a suitable Perzvat- 
sum and should avoid excessive distortion of the obli- 
que rays. Especially in the case of an objective of 
high relative aperture, aberrations calculated by 
Seidel theory will not be true aberrations, but the 
trend of their variations will, within not too wide a 
range, run approximately parallel with the true varia- 
tions, so that a few trigonometrical calculations will 
give an approximately linear relationship between the 
differences of the true aberrations and those determi- 
ned by SerpeL theory [3]. 

The following example of a triplet whose first lens 
is divided to lower the aberrations of the first surfaces 
will illustrate the method. To keep spherical aberra- 
tion low, the second lens is an Amici-meniscus and to 
ensure chromatic correction the first and last compo- 
nents are made of two glasses having the same refrac- 
ting index for yellow light, but different Anpe-num- 


bers. A high refractive index for the positive compo- 
nents and a low one for the negative gives a low PEeTz- 
VAL-sum and reduces the spherical aberration of the 
first. surfaces. 

The ratio of the distances of the positive thin lenses 
from the negative thin lens is taken as 2: 1, the first 
distance, for unit focal length, being d, = 0.5. If the 
contraction of the incident rays, w, for the negative 
lens, be 0.5, the power of the first thin lens, 9,, be 


; ; d; : : : 
unity, and taking 9,/9; = — to give low distortion for 


dy 
a diaphragm near the negative lens, it follows that 
3 = 2, 9, = — 2. The Petzvat-sum is then 0.355. 
The Amici-meniscus will have 9, = 0.494 and the 


first part of the positive system will have 9, = 0.575, 
if the separation between these thin lenses is 0.242. 
Using the convenient coefficients of ARGENTIERI, 
we get for the SempELt-sums YA, 2B, LP as used by 
Berek and for S,,0 the following formulae (by slide 
rule) : 
SA = 0.3055 wu; — 0.165 wu, — 
0.0739 p= 0125 use 
+ 0.0671 uf — 0.125 p, — 0.289 
> B = 0.452 ean = 
— 0.0827 us — 0.281 pp, + 
+ 0.1422 vs + 0.129 us — 0.0915 
XP = 0.0924 u{— 0.7867 up. + 
+ 0.3015 uy + 1.1035 us + 1.609 


where p = | : : F 
if if y 
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Four trials lead to SA 0.64, YB 0.32, Sr unre 
from the above formulae, which give for the SripEL- 
sums in areal system of lenses with suitable thicknesses, 
S045, 21, (S10) 0.92, SII-—0.074, STV 0/233: 
(The Serer calculations for the thick lenses require 
rather higher accuracy than is given by the ordinary 
slide rule). 

From this stage corrections are made by variations 
of the system of thick lenses as mentioned above. 
Six trials lead to a system of the following characte- 
ristics. 


Triplet H9m Pier == 
F. e! n’ Vv 
0.6671 0.2 1.7426 + 5.0 & 27.7 
12201 0.003 1 
0.5480 Ost 1.7425 + 5.0 
0.7102 0.1945 a 
—1.0702 0.03 1.4645 65.7 
0.3001 0.06 1 
0.4968 0.25 1.7426 AHO) & Die 
—1.0J353 0.3870 1 
Table of SEIDEL coefficients of the surfaces 
A B ify E P 
0.82381 0.54955 0.36659 0.67069 0.63882 
0.00326 0.02885 0.25541 —0.82815 —0.34897 
—0.00213 —0,00272 —0.00348 0.98924 0.77764 
0 0.00055 0.59915 —0.97104 —0.60004 
—0(0.75679 0.36813 —0.17908 0.23128 —0.29638 
—0.52440 —1.07038 —2.18477 —6.61643 —1.05679 
0.13029 0.39668 1.20770 6.28858 0.85782 
0.59658 0.00338 0.00002 0.00220 0.38907 
Sums 0.27061 0.27404 = 0.06155 —0.23363 0.36117 
Sie 8 027) > 2, (y T= 0) = 1.0126 
Se Oey eye = IAD Vie (0223 


Z, is chosen greater than before to get a better cor- 
rection of distortion. The coefficients in the first three 
columns of the upper table are low, ensuring a stable 
correction of spherical aberration, coma and curvature 
of field. The final Petzval sum is of the same order as 
before. Only the fourth column shows large coeffi- 
cients for the sixth and seventh surfaces, which 
means that distortion will not be so stable, but in the 
case of this rather asymmetric system which allows 
only limited field angles (up to + 18°), this is not of 
great importance. Figure 1 shows the final form of 
lens. 
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In figure 2, (a) shows the spherical aberration ; 
(b) shows the wave surface using coordinates such that 
spheres with their centres on the axis of the system 
are straight lines ; in this graph the abscissae are then 
nearly proportional to the areas of the corresponding 
zones of aperture ; and (c) is a diagram of the cosine 
of the differences in wavelengths between the actual 
wavefront and the closest fitting sphere. Figures 2 (0d) 
and (c) are drawn for a focal length of 25 mm [4], 
[5]. The last curve shows that the average cosine of 
the phase differences of the wavefronts is 0.7 so that 
the intensity of light at the focus will be 50°% of the 
intensity which could be reached by a perfect sys- 
tem ; this is considered tolerable. 
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The stabilised retinal image 


R. W. Ditcupurn and D. H. FEenpER 
University of Reading, England 


Summary. — A beam of light reflected from a mirror attached to a contact lens, worn by the subject, produces a lest object whose image 
does nol move across the retina in response to eye-movements. This is called the stabilised retinal tmage. In this paper, conditions 
for accurate stabilisation are discussed. It is shown that a small defect of stabilisation is important and that normal vision may be 
restored by introducing small movements which simulate some of the normal movements of a retinal image. Experiments on the 
stabilised image with light interrupled at different frequencies are described and it is shown that normal vision is obtained near 


the flicker fusion frequency. 


Sommaire, — Un rayon lumineur, réfléchi par un miroir fixé sur un verre de contac! porté par le sujet, forme Vimage @un lest objet dont 
Vimage est immobile sur la rétine quels que soient les mouvements du globe oculaire ; cette image est dite « stabilisée », Le présent 
article traile des conditions requises pour une stabilisation tres précise. Il sera montré qwun léger défaut de stabilisation est impor- 
tant pour la vision normale et que celle dernisre peut étre rétablie a aide de petits déplacements de Vobjet remplacant quelques- 


uns des mouvements normauc dune image rétinienne. 


Des expériences sur Vimage stabilisée, avec illumination intermiltente a différentes fréquences, sont décrites ; il est démontré 
que la vision normale se trouve rétablie pour une fréquence proche de la fréquence crilique de papillotement. 


ZUSAMMENFASSUNG. — Man kann ein sogenannles stabilisiertes Nelzhautbild in folgender Weise erzielen. Die Versuchsperson trdgt 
ein Kontaktglas, an dem ein kleiner Spiegel befestigt ist. Ein Lichtbtindel, das an diesem Spiegel reflekliert wird, dient zur 
Erzeugung des Bildes einer Pritfmarke, das trotz der Augenbewegungen auf der Netzhaut nicht auswandert. Die Bedingungen 
fiir eine genaue Stabilisierung werden untersucht. Es zeigt sich, dass chon ein kleiner Fehler in der Stabilisierung von Bedeutung 
ist und dass das normale Sehen durch kiinstliches Hervorrufen kleiner Bewegungen, welche die natiirlichen Augenbewegungen 
anndhernd nachahmen, wiederhergestellt werden kann, Es werden Versuche an dem stabilisierten Bild beschrieben, wo bei die Beleuch- 
tung mit wechselnder Frequenz unterbrochen wird. Es ergibt sich, dass in der Nahe der Flimmerfrequenz normales Sehen eintrilt. 


1. Introduction. — Earlier work on eye movements 
has shown that, when a subject is fixating as stea- 
dily as possible, the following eye movements are 
present 


(1) a tremor of amplitude of the order of 15 sec. 
arc and frequency ranging from 30 to 80 ¢. p. s. 

(ii) a series of “ flicks ” of up to 20 min. are occur- 
ring at irregular intervals ranging from 0.03 to 5.0 sec. 

(i) slow drifts in the intervals between flicks. 

Recent researches have shown that these eye move- 
ments play an important part in the visual process. 
Attempts have therefore been made to assess their 
importance by studying the effect on visual acuity 
when the eye views an object which moves so that 
its image on the retina does not move across the retinal 
receptors. This is described as a stabilised retinal 
image. 

The production of a stabilised retinal image was 
first reported by DircupurN and Ginspora [2] who 
noted that if the image consisted of a circular field 
divided vertically into halves of unequal brightness, 
the line of demarcation gradually disappeared. Shor- 
tly afterwards Rariirr [9] also reported the pro- 
duction of a stabilised retinal image ; subsequently 
Rices, Ratiirr et al. [10] gave the results of an exten- 
ded investigation in which they reported that the 
ability to see a fine black line on a bright field gra- 
dually diminished whilst viewing such an image. 


2. Methods of producing a Stabilised Retinal Image. — 
Dircupurn and Grinspore ([2] and [3]) used an 
apparatus as shown in figure 1. A test object is illumi- 
nated by parallel light. The lens L focuses the light 
on a mirror M worked on the side of the corneal bulge 
of a contact lens and forms an image of the test object 


on a screen. The eye views the image on the screen 
in the two mirrors M, and M;. 

These mirrors are arranged to give a visual path 
twice as long as the distance from the mirror M to the 
screen, so that the angular motion of the image seen by 
the eye is the same as that of the eyeball. The image 
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Fic. 1. — Original optical system for producing a stabilised 
retinal image as used by Dircupurn & Grnspore | 2 Ik 
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Fic. 2. — Improved optical system for producing an accurately stabilised retinal image. y M, 


on the retina therefore remains nearly stationary 
when the eye rotates in a horizontal plane. Stabilisa- 
tion may be regarded as complete if the residual mo- 
tion of the image on the retina is less than the size 
of one cone. This corresponds to about 0.1 min. are 
in the visual field. With the apparatus just described 
stabilisation would not be complete in the vertical 
plane even if it was complete in a horizontal plane. 
Since the mirror M must be worked on the side of the 
corneal bulge of the lens (the normal to the mirror 
makes an angle of 35° with the visual axis in the lens 
used by Ditrcupurn and GinsBore [2]) only a compo- 
nent of the upward rotation of the eyeball appears in 
the upward movement of the reflected beam. With the 
system described above, the angular motion of the 
image in the vertical plane is only about 85 % of the 
motion of the eyeball. This fact led Dircupurn and 
Ginspora [2] and Riacs and Raruirr [10] to use 
test objects with vertical patterns, so minimising the 
effect of lack of vertical stabilisation. 

Another factor tends to reduce the stabilisation in 
both planes. It has been shown by Fenper [5] that 
the eye makes torsional motions (that is, rotation 
about the visual axis) when fixating. These motions 
are of a random nature, not intimately connected 
with the other motions, and may range up to about 
10 min. arc. If the normal to the mirror M is not paral- 
lel to the visual axis, components of this torsional 
motion appear in the vertical and horizontal motion 
of the projected image. These components may be as 
large as 5 min. arc, that is, sufficiently large to des- 
troy the stabilisation completely, 


If the effects mentioned above are to be avoided, the 
mirror M must be positioned so that the normal to 
its plane is parallel to the visual axis. This is done by 
mounting the mirror on a stalk attached to the con- 
tact lens as described by Fenner [5]. The mirror is 
fixed to a small steel ball which fits into a socket at 
the end of the stalk. This allows the mirror to be 
fitted to the contact lens after the latter has been 
placed in the eye. The plane of the mirror can then be 
accurately adjusted. The use of the stalk has the addi- 
tional advantage that the mirror is not dimmed by 
eye fluids. When the subject is wearing a contact 
lens and the other eye is covered, the tendency to 
blink is much less than usual. Trouble due to distur- 
bance of the system by blinks is very much less than 
might be expected and does not constitute a serious 
difficulty. 

One further effect which impairs stabilisation is due 
to the motion of the mirror itself. Rotation of the 
eyeball causes both rotation and translation of the 
mirror attached to the contact lens. If the transla- 
tory motion has a component in the direction of the 
incident beam of light, it will cause a corresponding 
shift in the projected image. In the system shown in 
figure 1 this may amount to about 1 min. arc, which is 
again enough to destroy stabilisation. To avoid this, 
the projected beam of light should be incident on the 
mirror along the line passing through the centre of 
the mirror and the centre of rotation of the eye. 

An apparatus designed with regard to the preceding 
discussion is shown in figure 2. Light from the ball of 
a pointolite lamp is rendered parallel by L, and illu- 
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minates the test object. L, focuses the source on the 
mirror M, and projects an image of the test object 
at lie 

The whole of this projector is for convenience pla- 
ced above the level of the subject’s eye; the light is 
brought down to the correct level by a periscope 
arrangement, the lower reflector of which is turned 
so that light strikes M, in the direction recommended 
above. 

L, forms an image of I, on the ground glass screen 
S by reflection at M, and Mg. L; reverses the direction 
of motion of the image with respect to the eye but 
this is corrected horizontally by the pair of mirrors M, 
and M, and vertically by the Dove prism. 

In these circumstances, the retinal image is sta- 
bilised when. 

(1) 
where 

g = equivalent path length from L, to 5, 

f = focal length of Lg. 

% = distance from mirror M, to first image I,, 

| = distance from nodal point N of the eye to the 
screen S. 


2(qq—fz=lf 


Fic. 3. — Effect of nodal point and centre of rotation of the 
eye on Stabilisation. 
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The use of an image transmitted through a ground 
glass screen enables much higher brightnesses to be 
obtained than when the image is seen by reflection 
from a screen. For work at even higher brightness 
levels, the screen S may be replaced by a lens of focal 
length SI, (I, is the second image of the pointolite 
ball), the image is then seen by Maxwellian view. 

The position of the nodal point and the centre of 
rotation of the eye must be known in order to set up 
the stabilising system accurately. This is shown by 
the following calculation. 

Referring to figure 3, if the visual axis turns through 
an angle 6 from the position VA to VA‘ about the 
centre of rotation O, then the projected image moves 
from T to T, where TT, = 2 6d. d is the distance 
from the mirror to the image. This image is repeated 
in the plane SS’ either by the mirror system of figure 4 
or the optical system of figure 2. In either case, the 
effect is that the image moves from T’ to 1, (equal 
to TT,) SS’ being at a distance D from the nodal 
point of the eye; the magnitude of D is determined 
below. If the nodal point N is at a distance r from 
the centre of rotation, then N will move to N’ where 
NN’ =r 0.The linear motion of T’ relative to N is 
thus (2 d 6 —r 0) and the relative angular motion 
is 0 (2 d —r)/D. Call this angle o. 

Let R be the distance from the retina to the centre 
of rotation of the eye; the nodal point is thus at a 
distance (R + r) from the retina. The linear motion 
of the retinal image relative to N is (R + r) @ or 
(R +r)s—r0do relative to the original image 
point P. However this point P on the retina moves 
to P’ where PP’ = R 0, thus the motion ¢ of the 
image relative to the retina is given by 


— } (R+r)o—ro{—Ro 
(RB 7) (38) 


If the retinal image is to be stabilised, « must be 
zero, hence (g — ) must be zero. 


But ( pe care) 
D 
_ 6(2d—r—D) 
D 
Thus (2 d —r — D) = 0 
or 
(2) D = (2d—r) 


Unfortunately the positions of the nodal point and 
centre of rotation when wearing a contact lens are 
not easily determined. However, the optical system 
of figure 3 may be set up assuming likely positions 
for these two points and then a fine adjustment may 
be made in the following way. 

The mirrors M, and Mg are carried on a stage sli- 
ding along an optical bench parallel to the visual axis. 
_ The subject is asked to view a test object consisting 
of a black line crossing a bright circular field. Initially 
the retinal image is only partially stabilised, and al- 


vol 2;.n° 8, Oct. 1955) STABILISED 


Z Nasalward 


Motion of Image 


Fic. 4. — Effect of movement of the nodal 
point on stabilisation. 


though the line disappears, it quickly reappears, pre- 
sumably because motions of the eyeball can shift the 
image to new parts of the retina. The subject presses 
a key when he can see the line ; this operates a timing 
mechanism and it is usually found that the line is 
visible for about 95 °% of the viewing time. The stage 
carrying the mirrors is then moved slowly along the 
optical bench ; this varies the factor g in equation (1) 
and so alters the magnification of the optical system, 
which has the same effect as varying D in equation (2). 

At one position, it is found that the percentage of 
time for which the line is visible drops dramatically, 
being as low as 50 % for a black line of width 5 min. 
are on a field of brightness 25 millilamberts. In the 
apparatus in use at present, it is found that a varia- 
tion of g by about + 2 mm. in 66 cm. reduces this 
stabilisation effect considerably. It is possible that the 
precision with which g can be determined by this 
method will lead to a measurement of the geometry of 
the optical system of the eye. 

The nodal point of the eye does not occupy a 
unique position. It is known to change with accommo- 
dation and it has been suggested by FENDER [6] that 
the position of the nodal point is sensitive to colour. 
This, combined with the fact that the optical system 
of the eye is used off-axis, means that a change in 
accommodation can move the image across the retina 
(see fig. 4) and so destroy stabilisation. 

The effect of this movement may be reduced by 
arranging that the image movement on the retina 
causes important straight boundaries in the pattern 
to move along their length. The angle « (fig. 4) is mea- 
sured for each subject and the corresponding angle 
in the vertical plane. The boundaries in the pattern 
are then set to make an angle tan (8/«) with the 
horizontal. 

Arnutr [1] and co-workers have noted periodic 
fluctuations of accommodation, which occur at 
about the same frequency as the disappearance and 
regeneration of a stabilised image. These fluctuations 
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however are so small (of the order of 0.02 diopter) 
that they alone are probably not responsible for the 
regeneration of a stabilised image. A more probable 
explanation is that the subject, being presented with 
a most unusual visual experience, makes some convul- 
sive movement and so upsets the stabilisation of the 
image. It is noticeable that older subjects with small 
accommodative range experience more complete sta- 
bilisation than do young subjects. 


3. Effects of a stabilised Retinal Image on Vision. — 
The most remarkable effect of stabilisation is of 
course the very rapid disappearance of a_ stabilised 
retinal image. This has been dealt with very fully by 
the authors mentioned, but one point should be em- 
phasised — this is not a fatigue effect in the normally 
accepted sense; for if a black bar crossing a bright 
field forms the stabilised image, then, as reported by 
DircuBuRN [4] it is the black bar which disappears ini- 
tially, the whole field becoming of uniform brightness. 

The manner in which large illuminated patches 
disappear is of interest ; circular patches first become 
diffuse at the edges, the effect then spreads both 
inwards and outwards, the image gradually becoming 
less bright and finally disappearing. Sometimes this 
effect is not complete and the whole visual field re- 
mains suffused with light. If the test patch has sharp 
corners, these are always the first to fade, reducing 
the image to a roughly circular patch which then 
behaves as described above. 


4. Effects of introducing Artificial Movements to a 
Stabilised Retinal Image. — The function of eye 
movements has been studied by causing a stabilised 
retinal image to move in a fashion which simulates 
the effect of eye movements ; in every case it is found 
that an immediate regeneration of the image occurs. 

Firstly, it has been found that tapping the appara- 
tus introduces small vibrations which move the stabi- 
lised image through about 1 min. arc at a fairly high 
frequency, possibly about 100 c. p. s. such a vibration 
is similar to the ‘‘ tremor ” of a normal eye move- 
ment. This immediately regenerates a stabilised image 
if it has disappeared and also assists the eye in focus- 
sing on the test image. It appears that the information 
which enables an eye to focus an image is lost when 
the image is stabilised. More precise experiments of 
this kind are in hand at present. 

Secondly, it is arranged that when a_ stabilised 
image (consisting of a fine dark line on a bright field) 
disappears, the subject presses a key. This immedia- 
tely moves the image by about 10 min. arc on the 
retina, thus simulating a flick. Again the image is 
immediately regenerated It is not suggested howe- 
ver that this is the primary function of a flick. In this 
respect the explanation offered by DircuBurn and 
GinsporG [3] seems much more probable. 

5. Experiments on a Stabilised Image with inter- 
rupted Light. — In previous experiments with the 
stabilised image, the illumination was continuous. 
We shall now describe some experiments with inter- 
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I'ic. 5, — Visibility of a stabilised retinal image 
in interrupted light. 
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rupted light. In one experiment the subject viewed a 
test object consisting of a black line 5 min. are wide 
on a bright circular field, diameter 1° and brightness 
25 millilamberts. The source of light was then inter- 
rupted (equal periods of dark and light) at various 
frequencies. The subject pressed a key which opera- 
ted a timing device whenever the black line was 
visible. Figure 5 shows the percentage of time for 
which the line was visible at several frequencies. The 
total viewing time was always 30 sec. It will be noti- 
ced that the curve has a peak of 97 % at about the 
flicker fusion frequency for the subject. Points on 
the curve could be determined with a probable error 
of about + 5% except the peak value, which is repro- 
ducible to a higher order of accuracy. 
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Whilst conducting these experiments, the subjects 
complained of the brightness of illumination at high 
interruption rates. This was surprising as the mean 
illumination remained constant at about 10 milli- 
lamberts throughout the experiment. One subject 
in fact estimated the brightness to be of the order of 
1000 milllamberts ! 

In order to test this effect, an attempt has been 
made to measure the apparent brightness of the inter- 
rupted image. The projector of figure 3 is replaced by 
that shown in figure 6. The stabilised image now con- 
sists of a circular field divided into halves. One half is 
illuminated with continuous light and the other with 
interrupted light. A photometric wedge controlled 
by the subject is introduced into the interrupted 
light path and the subject is asked to match the bright- 
ness of the two halves of the field whenever the conti- 
nuous half regenerates. 

Some very anomalous results have been obtained. 
Initially the interrupted half of the field is set by the 
experimenter to a lower brightness than the conti- 
nuous half; the subject describes the interrupted 
half as the brighter, but when told to make a match, 
always moves the wedge so as to make this half brigh- 
ter still. This of course is the direction in which the 
wedge should be moved to achieve a true photometric 
match between non-stabilised images, but the sub- 
ject always overshoots the correct position and con- 
tinues searching to the end of the wedge. 


6. Discussion. — Work on the electrical impulses 
transmitted by individual fibres of the optic nerve 
(HARTLINE [8] and Granit [7]|) shows that, for some 
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retinas, a nerve fibre carries a large burst of impulses 
at a high initial rate when the receptor connected to 
the nerve is illuminated. If the illuminations is main- 
tained at a steady value, the rate of transmission of 
impulses falls to a low value sometimes to zero. 
There may be a further large burst of impulses when 
the illumination ceases. The exact behaviour of the 
individual receptors of the human retina is not known 
but it is likely to be generally similar to that descri- 
bed above. Dircusurn [4] has shown that this leads 
to a simple theory of the function of normal eye- 
movements and of previously reported results on the 
stabilised retinal image. 

The new observations described in this paper appear 
to fit in with this hypothesis. Let us consider a simple 
pattern consisting of a light area and a dark area. 
In normal vision the nerve fibres connected to recep- 
tors near the retinal image of the dividing line are 
strongly stimulated, because these receptors are 
subject to rapid changes of illumination owing to 
the tremor motion of the eyeball. The other move- 
ments of the eye ensure that no one set of receptors 
has this strong stimulation for very long because the 
image of the boundary is continually being moved. 

With continuous illumination and an accurately 
stabilised retinal image, all receptors receive constant 
stimulation and the image fades as the initial burst 
of nerve impulses dies down. The subject can no lon- 
ger appreciate a sharp edge and therefore it becomes 
difficult to focus. 

If stabilisation is not very accurate, the image is 
quickly regenerated either by a large fluctuation of 
accommodation or by the flick movement. The size 
of a flick ranges up to 20 min. arc, so that if the com- 
pensation of eye movements is only 99 % complete, 
a flick could move the image by about one cone dia- 
meter. Thus the adjustment of the apparatus for com- 
plete stabilisation should be fairly critical and the per- 
centage of time for which the image is not seen should 
increase rapidly as we approach complete stabilisa- 
tion. 

The extent to which the stabilisation described in 
this paper is better than that previously obtained is 
shown by the fact that Riacs and Raruirr [10] 
observed the disappearance of a line 35 sec. are wide 
on a field of brightness 5 ft-lamberts. In the present 

experiments a test object consisting on a line 5 min. 
arc wide on a field of brightness 25 millilamberts 
disappeared for 50 % of the time. 

With interrupted illumination and stabilisation, the 
same group of receptors is continually subject to chan- 
ging illumination. These receptors may thus, with 
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only moderate illumination, be made to generate 
repeated large bursts of impulses in the nerve fibres. 
This gives the impression of painfully bright illumi- 
nation. It may be that, for some nerve fibres, a satu- 
ration level is reached under these conditions. When 
this happens the subject, using the apparatus shown 
in figure 6, may be expected not to be able to judge 
whether the movement of the photometric wedge 
increases or decreases the illumination. 

The results shown in figure 5 can be explained in the 
following way. At very low interruption frequencies, 
normal vision is restored momentarily each time the 
illumination increases and probably each time it goes 
off. This should reduce the percentage of time for 
which the stabilised image disappears. The condition 
improves as the interruption rate increases until the 
time between successive periods of illumination is too 
small for the receptor/nerve fibre system to recover 
appreciably, whereupon the condition will worsen 
once more. When the frequency of interruption is very 
high, the illumination will be effectively continuous. 
From the above discussion, we expect the maximum 
of the curve of figure 5 to occur at about the flicker 
fusion frequency. We also expect the light to be effec- 
tively continuous when the frequency of interruption 
is of the same order as the rate of the initial discharges 
of the nerve fibres, that is, rather more than 50 c. p. s. 

It is, of course, likely that the above simple discus- 
sion will need a great deal of extension and modifica- 
tion when more detailed results are available. 
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Peripheral and foveal dark-adaptation curves as a function 
of luminance and wavelength (') 
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SUMMARY. 


Peripheral and foveal dark-adaptation curves were measured after stimulation by luminances ranging from 1300 fo 


0.7 cd/m? for red, yellow, green and violet light-adaptation field and test spot. It was found that the dark-adaptation time of the 
cone apparatus depends but slightly on the colour of the stimulating field and the test spot. The dark-adaptation time of the rod 
apparatus appeared to be strongly dependent on the colour of the stimulating field (fig. 4). Typical peripheral adaptation curves 
are shown in figures 1, 2 and 3. Some tupical phenomena in ‘oreal adaptation are mentioned in section Sips 


Sommatre, — Les courbes d’adaptation & Vobscurilé pour la fovéa et la région extrafovéale ont été délerminées aprés stimulation par 
des luminances allant de 1 300 & 0,7 cd/m? pour les lumiéres rouge, jaune, verte ef violelte du champ d’adaptation et de la plage 


de mesure. 


On a trouvé que le temps d’adaptation a Vobscurité des cones ne dépend que peu de la couleur du champ d’ adaptation et de la plage 
de mesure. Le temps d@’adaptation des bdtonnets dépend au contraire fortement de la couleur du champ d’adaptation (fig. 4). Des 
courbes types d’adaptation périphérique sont données sur les figures 1, 2, 3. Quelques phénoménes typiques obtenus dans l’adapta- 


lion fovéale sont mentionnés dans la section 3.2. 


ZUSAMMENFASSUNG, — Es wurden die peripheren und fovealen Dunkeladaptationskurven nach der Reizung mit Leuchldichten zwis- 
chen 1 300 und 0,7 cd'm? fiir rotes, gelbes, griines und violettes Lichtim Adaptationsfeld und Priiffleck gemessen, Es ergibt sich, dass 
die Zeit der Dunkeladaptation des Zapfenapparates nur wenig von der Farbe der Helladaptation und des Priffleckes abhdngt. 
Die Zeit fiir die Dunkeladaptation des Stébchenapparates hdngt aber stark von der Farbe des Helladaptationsfeldes ab (Abb. 4) 
Die Abbildungen 1, 2 und 3 zeigen typische periphere Adaptationskurven. Einige auffallende Erscheinungen bei der fovealen 


Adaptation werden im A bschnitt 3.2 erdrtert. 


1. Introduction. — The process of dark adaptation 
has been frequently investigated by measuring dark- 
adaptation curves, which represent the minimum per- 
ceptible luminance of a test spot against a dark back- 
ground as a function of time after exposure of the eye 
to a bright stimulating light source. The influence of 
the wavelength of the stimulating light upon peri- 
pheral dark adaptation was measured by HuxzBurt [1], 
the influence of the wavelength of the test light by 
KouLrauscH [2]. MANDELBAUM and Mintz [3] mea- 
sured the influence of both factors upon the foveal 
adaptation curve. 

The present article deals with an extension of these 
measurements. Peripheral dark-adaptation curves were 
measured 7° temporally from the fovea of the left 
eye of two observers (J. M. and H. H.) with mono- 
chromatic violet (486 my), green (546 muy), yellow 
(578 my) and red (644 my) stimulating light and mono- 
chromatic violet, green and yellow test light. In addi- 
tion, a non-monochromatic red test light with a domi- 
nant wavelength of 650 my was used. The monochro- 
matic light was produced by a high-pressure mercury- 
cadmium vapour lamp, combined with colour filters, 
the non-monochromatic light by a tungsten ribbon 
lamp with a temperature of 2 000°K, combined with a 
Corning glass filter C 2403. The luminance of the sti- 
mulating field was nearly 1000, 100, 10 and 1 cd/m? 
in different measuring cycles. Foveal dark-adaptation 
curves were measured for the same monochromatic 
violet, green and yellow stimulating and test light. 


(1) The investigations described in this article were made 
under the auspices of the ‘‘ Nederlandse Stichting voor Ver- 
lichtingskunde ”’. The author is now with Research Labora- 
tory AKU and Affiliated Companies, Arnhem. 


2. Peripheral measurements. 


2.1. Measuring arrangement and procedure. — 
The measurements of peripheral dark-adaptation 
(d.a.) curves were made with the aid of an adaptome- 
ter, consisting principally of a circular light-adapta- 
tion (I. a.) field with a diameter of 14° and a square 
test spot with sides of 2°. The latter was situated in a 
direction perpendicular to that of the |. a. field, but 
could be seen by the observer in the direction of the 
l. a. field with the aid of a rotatable mirror. The spot 
of the retina on which the 1. a. field and the d. a. test 
spot were focused was fixed with the aid of a small red 
fixation light. The luminance of this fixation light 
was dimmed gradually during dark adaptation. The 
observer used an artificial pupil with a diameter of 
2.4mm ; his head was held in one position by a mouth- 
piece mounting. The test spot was shown to the obser- 
ver during 0.18 sec each second. 

The luminance of the |. a. field was adjusted by 
means of grey filters to a number of values, ranging 
irom 1300 to 0.7 cd/m?. The luminance of the d.a. 
test spot could be continously adjusted by the obser- 
ver during dark adaptation by moving a diaphragm 
in the light beam which illuminated the test spot. The 
setting of the diaphragm was recorded on a strip of 
paper as a function of time. When the luminance of 
the test spot reached the instantaneous brightness 
threshold value, the observer pushed a button, making 
a mark on the recording strip. 

After a pre-dark adaptation of half an hour, the 
observer adapted to the l. a. field for 3 minutes. Then 
he generally measured successively three adaptation 
curves for different |. a. luminances or colours. The 
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first measurement of the luminance threshold after 
I. a. occurred 12 seconds approx. after the beginning 
of the d. a. Afterwards, as many threshold values were 
determined as the observer thought necessary for a 
correct determination of the d. a. curve. In the period 
of quickly changing adaptation conditions, the thres- 
hold was considered to be that luminance of the test 
spot at which this spot became just invisible for the 
observer. In the period of slowly changing adaptation 
conditions, the luminance at which the observer esti- 
mated the chance of observation of the (periodically 
shown) test spot to be 50% was taken as threshold 
value. The luminances of the stimulating field and the 
test spot were measured with the colour-filter me- 
thod ([4] p. 297), using a photometer with Maxwellian 
view ({4] p. 152). One half of the lens of the photometer 
projected light of the stimulating field or test spot to 
the observer’s eye, the other half light of a standard 


cd/m? 


Fic. 1. — Peripheral dark-adaptation curves measured by 
J. M. with a violet monochromatic test light (Aq = 436 my) 


after stimulation by monochromatic violet light (A, = 436 my) 
for 4 luminances B,. By, is the minimum brightness thres- 


hold of the cone apparatus, B,,, that of the rod apparatus. 
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Fic. 2. — Peripheral dark-adaptation curves measured by 


J. M. with a monochromatic violet test light (Aq = 436 my) 


after stimulation by nearly equal luminances B, in cd/m? 
and monochromatic light of wavelength A, (in my). 


source (tungsten ribbon lamp of known tempera- 
ture). 

The measurements with green and yellow light 
were checked with a calibrated Weston photronic cell 
with Viscor correction filter. The agreement between 
the two methods was better than 10%. 


2.2. Measuring results. — Only the results of obser- 
ver J. M. will be given, as those of H. H. are almost 
identical. 

Adaptation curves showing the general aspects of a 
change of the luminance of the stimulating field, the 
colour of that field and the colour of the test field are 
presented in figures 1, 2 and 3 respectively. The break 
in some of the curves indicate the presence of both 
cones and rods peripherally. The minimum luminance 
thresholds of the cone and rod systems are indicated 
in the figures by Bme and Bmr respectively. 

The curves of figure 2 were shifted along the 
ordinate axis so that the thresholds Bm, coincide. 
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Fic. 3. — Peripheral dark-adaptaticn curves measured by 
J. M. after stimulation by monochromatic violet light 
(Ay = 436 my) and the luminance B, = 95 cd/m? with four 
test spots, emitting light with wavelength Aq. 


The actual values of the thresholds for wavelength 
436, 546, 578 and 644 mu were 1.4, 0.85, 1.41 and 
1.2 x 10-5 cd/m? respectively. 

The different values of B,, for different A, in 
figure 3 are due to the unit of luminance used; viz. 
the cd/m?. This unit is related to cone vision. However, 
the second part of the curves of figure 3 refers to rod 
vision. If the ordinates were expressed in ‘‘ scotopic ” 
cd/m?, the values Bes calculated from all measure- 
ments with the aid of the scotopic luminosity factor 
V5, (CIE, 1951) would differ by less than a factor of 
2.5 (table I). The formula used in this calculation can 
be written : 

p46 LNs, 

(1) Bose a 680 Vy Brar 
where the luminances are assumed to be monochro- 
matic. 


Adaptation times t;,, were derived from the adap- 


tation curves measured. ¢),, is the time required to 
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observe a luminance Bi, which is 1/4 log unit (a fac- 


tor 1.8) greater than the minimum luminance thres- 


TABLE I 


Average threshold luminance, expressed in photopic cd/m? 
(By) and ‘scotopic’ cd/m? (B,,,), derived from all measu- 
rements and the CIs Visibility curves for photopic and 
scotopic vision (Vj, V3 resp.). 


rt if 
A in mu [Bie V/V. [Biro 
| 
650 3.8 x 102 0.00632 6.2 x 104 | 
578 DB se AOE 0.153 9.0 x 10-4 
546 yy Se i=“ 0.557 3.6 xX 10-4 
436 | 122 052 15.3 4.6 x 10-4 


hold Bmr. This time was considered to be the maxi- 
mum adaptation time which could be derived with 
sufficient accuracy from the adaptation curves. The 


are collected in table II. 
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TABLE II 


Dark adaptation time Ola for the rod system as a function of 


the wavelength A, and the luminance B, of the stimulating 
field and the wavelength Ag of the dark-adaptation test 
light. 


Ag in mit 
ne Byin cd/m? 436 546 578 650 
in mit a0: a eae 
tyyin see uygin see oe in sec vain sec 
se alee se Ore! elm Or 
436 13.0 a0 32 22 is 
0.95 15.0 10.5 9.5 Thee 
0.095 7.0 7.0 4.6 Bes 
0.0070 Sr 3.0 
546 T3E2 oo UBS) Or 7.8 
1e20 5.6 Haz 4.5 2.6 
0.108 2.6 Dee 1.65 0.85 
578 13.5 7.8 7.8 5.6 4.6 
usa es} AnD, 3.8 PsA 1.40 
0.130 1.85 Lao 1230) 0.85 
644 Tigh Qed 2.6 Dro 11 PAs) 
0.75 (1.4) 1.45 150) 0.90 
0.092 (122) 0.95 125 0.64 


Cc 


Similarly, adaptation times t,,, were determined 
from the adaptation curves in which a cone plateau 
could be determined. Here ¢;, is the time required 


to observe a luminance Bi), which is 1/2 log unit 


(a factor 3.2) greater than the minimum luminance 
threshold Bme of the cone apparatus. These data are 
collected in table IIT. 
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TABLE III 


Dark adaptation time i 2 for the cone system as a function 
of B,, A, and dg. 


Aq in mu 
MK Byin cd/m? 436 546 578 650 
in mu < 10-2 
€ ¢ Cc Cc 
1/2 te tye Mle 
in sec in sec in see in see 
436 13.0 110 66 46 85 
0.95 26 14.0 10.0 9.0 
0.095 (6.0) (3.0) 
546 a2: 21 20 24 18.0 
1.20 PA) 
| GWisjt|| ile yes 18.0 22 22 
| 1.18 6.0 
2.3. Discussion. — Table II shows that the adap- 


tation time of the peripheral retina is greatest for a 
violet test hight and decreases for a green, yellow and 
red test light respectively when 2, and B, are taken 
as constant. The qualitative course of the curves for 
hive as a function of 2, is almost the same for the diffe- 
rent wavelengths 4,, when B, is taken as constant. 


Therefore, the mean value ¢,),, of an for different 
wavelengths 2, at the same wavelength A, was calcu- 
lated and plotted against 4, for luminances B, = 1 000, 
100 and 10 cd/m? (fig. 4). 


400 500 600 ML 700 


Frc. 4. — Mean dark adaptation time /,), for the rod system 


after stimulation with monochromatic light of wavelength A, 
at a luminance B, of 10, 100 and 1 000 cd/m? resp. 


FOVEAL DARK-ADAPTATION SW, 
A similar curve to that found for B, = 10 ed/m? 
was measured by Hutsurtr [1]. Because of diffe- 


rences in measuring procedure, the curves are not 
fully comparable ; the slope measured by the present 
author is larger than measured by Hutsurt but in 
better agreement with measurements of Hecur and 
Hs1a [5]. The latter authors measured a factor 3 for 
the ratio of adaptation time after stimulation by 
white light and red light respectively. 

The curves in figure 4 were drawn for B, = cons- 
tant (10, 100 and 1 000 cd/m? resp.). Re-calculating the 
measured results in scotopic luminances with the aid 
of eq. (1), it appeared that the curve for B, = 100 sco- 


topic cd/m? (*) is nearly independant of the wave- 


length %,. The constant value of ies found in this way, 
is 420 sec approximately. Because of small influences 
of A, (table II) and the different shapes of the curves 
of figure 4 at the longer wavelengths, it is not possible 
to generalise the conclusion drawn from the curve 


for B, = 100 scotopic cd/m? to all other values of Be 
From table III it can be seen that the influence of 


green and yellow stimulation light upon ti is small. 
The same result was found by ManpreLBaum and 
Mintz [3] in foveal measurements. That, relatively, 
ae after stimulation by violet light is greater than 
for the other wavelengths and also greater than was 
found by ManpELBauM and Min‘z, can be understood 
from the fact that at 7° peripherally macular pigment 
is absent. However, M. and M. did not find such a 
great influence of the colour of the test spot upon 


ti, after stimulation by violet light as is shown in 
table II. 


3. Foveal measurements. 


3.1. Measuring arrangement and procedure. — 
Except for the following items, the measurements 
were made in a similar way to that described in 2.1. 

The sides of the square test field subtended a visual 
angle of 1°. The red fixation light was replaced by an 
arrangement shown in figure 5. The fixation point is 
the point of intersection of the two rows of fixation 
lights. During light adaptation the observer looked 
as far as possible at the centre of the circular stimu- 
lating field. 

The minimum luminance at which the test spot 
could be observed was determined in three different 
ways : 

a. similar to that described in 2.1, but with the 
test spot shown for 0.033 second every 2 seconds ; 

b. the test spot was shown continuously. The obser- 
ver determined the luminance at which the test spot 
just disappeared ; : 

c. the observer determined the chromatic threshold, 
that is the luminance of the test spot at which he 


I 
(1) This value was chosen because the values of B, then fall 
in the range of values given in fig. 4. 
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Fic. 5. — The fixation point F is defined as the imaginary point 
of intersection of the lines of point sources ab. The points 
b were presented at an angle to the fixation axis of about 
6° and the points @ at an angle of 13°. 


could just detect the colour of that spot. The colour 
was previously unknown to him as the test colours 
were changed in a random order and the luminance 
was only increased starting from an invisible test spot. 
The observed colour was checked afterwards. 

The pre-dark adaptation time was decreased to one 
quarter of an hour. Each dark adaptation curve was 
derived from two to four runs. 


3.2. Measuring results. — Most measurements were 
made with violet and yellow stimulation light and 
test spot, in order to make conditions for adaptation 


Cc 
h 


as extreme as possible. The adaptation time t),,, 


defined in 2.2, was measured for continuous test spot 
and for chromatic thresholds (Cf. 3.1 b,c). The dark 


adaptation time ¢,), after stimulation by the same 


colours and aluminance of 1 300 cd/m? differed only 
by a factor 1.5 for the different wavelengths of the 
test light. The significance of the differences was not 
shown, because the measurements were difficult and 
rather inaccurate. The average time 1), after stimu- 
lation by violet light was 1.3 times greater than after 
stimulation by yellow light. Generally, the results did 
not contradict those of MANDELBAUM and MinvTz. 

It was found that the foveal dark-adaptation curve 
measured by J. M. after stimulation by yellow light 
(B, is 1.9 x 104 cd/m?) had a break both for periodi- 
cally and for continuously presented violet test spots 
(fig. 6). The break might be caused by the presence of 
rods in the fovea (PoLyak, [6]). However, there is a 
difference with regard to the peripheral results, be- 
cause no such a break was found when instead of 
yellow light, green stimulation light was used with a 
stimulation luminance up to 13.104 cd/m? (fig. 6). 
The possible influence of rods in foveal adaptation 
measurements complicates an interpretation of re- 
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Fic. 6. — Foveal dark-adaptation curve measured by J. M. 


after stimulation by 

1. monochromatic yellow light (A, = 578 mu) and lumi- 
nance B, = 1.9 x 104 cd/m? with a periodically presented 
violet test spot (Aq = 436 mu), subtending a visual angle 
of 1° approx. The measured points were obtained at 3 mea- 
suring sessions which are shown differently. 

2. monochromatic green light (Ay = 546 my) and lumi- 
nance B, = 13.104 cd/m? (—.—). Average curve from 


3 measuring sessions. 


Fic. 7. — Foveal adaptation curves measured by J. M., after 
stimulation by monochromatic violet light Q, = 436 mu.) 
and luminance B, = 14 x 10? cd/m?, with a monochroma- 
tic violet test light (Aq = 436 mu) which was presented 


either continuously (circles) or periodically (0.033 sec. per 
2 sec., dots). 


ane curves were measured in two separate measuring 
cycles. 
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sults such as those of ManpELBAuM and Mintz upon 
the basis of a trichromatic colour system. 

The adaptation curves measured with a periodi- 
cally presented violet test spot (Cf. 3.1 a) after stimu- 
lation by violet light and a luminance B, of 14 x 102 
cd/m? appeared to be different from those measured 
after the same stimulation but with a continuously 
presented violet test spot. The ratio of the luminance 
threshold for a periodically presented test spot to 
that for a continuously presented test spot increased 
with dark adaptation time. This ratio increases from 
! to 4.3 for the curves shown in figure 7 and is of the 
order of magnitude to be expected from BLONDEL 
and Rey’s formula [7], extended for supra-liminal 
light flashes by Hampron [8]. 
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Sur la sensibilité différentielle de l’ceil aux couleurs dans l’échelle de Newton 
notamment pour les teintes sensibles 


Erik INGELSTAM ET Lars JOHANSSON 


Laboratoire d’Optique, Ecole Royale Polytechnique Supérieure, Stockholm 70 


SoMMAIRE. — Une étude expérimentale a été faile sur la sensibililé différentielle de Veil pour les couleurs interférométriques selon 


Véchelle de Newron, a Vintérieur et autour de la teinte sensible du premier ordre. Le but de cette étude était la mesure des gradients 
du chemin optique au moyen dune lame double de quartz envisagée par FRANGON, le champ de vision étant par conséguent en 
forme de ruban assez étroit. L’écart moyen des mises au point se révéla étre de 8 A environ dans le domaine le plus sensible, Un 
gradient ou une différence de chemin optique quelconque peut étre mesuré avec la précision correspondante de 20 A, au moyen d’un 
compensateur fonclionnant comme méthode de zéro 4 la teinte sensible. Les coefficients de diffusion propre des liquides trés dilués 
sont accessibles avec la méme précision que ceux des liquides plus concentrés mesurés avec les méthodes ordinaires. 


Die Fahigkeit des Auges, die interferometrischen Farben in der Newtonschen Skala zu unterscheiden, beson- 
ders an und neben der empfindlichen Farbe erster Ordnung, wird einem experimentellen Studium unterzogen. Es sollen dabet 
Gradienten der optischen Wege mit Hilfe der von FRANGON angegebenen Quarzdoppelplatte gemessen werden, In Ubereinstim- 
mung hiermit ist das Blickfeld bandférmig angeordnet. Die entsprechende mittlere A bweichung der Einstellungen wurde in dem 
empfindlichsten Gebiet zu 8 A bestimmt, Ein optischer Gradient oder eine Gangdifferenz irgendwelcher Grdsse kenn mit der 
entsprechenden Genauigk7it von 20 A gemessen werden, wenn man einen mit der empfindlichen Farbe als Nullmethode arbeiten- 
den Kompensator verwendet. Selbstdiffusionskoeffizienten von s°hr verditnnten Lésungen kénnen so mit derselben Genauigkeit 
gemessen werden, wie man sie bei konzentrierteren Losungen mit gewéhnlichen Methoden erhalt. 


SumMMARY. — An experimental determination was made of the eye’s ability to discriminate interference colours along the Newtonean 


scale within and on both sides of the first order sensitive colour. The application aimed at is to measure optical path gradients 
by means of the double quartz plate as suggested by FRANGON and accordingly the field of view consists of a fairly narrow band. 
The mean deviation of the settings was found to correspond to about 8 A in the most sensitive region. Any gradient or path diffe- 
rence value can be measured with the corresponding accuracy of 20 A by means of a compensator acting as zero method at the 
sensitive colour. Coefficients of self-diffusion of very dilute solutions may be measured as accurately as those of higher concentration 


with ordinary methods. 


Introduction. — En interférométrie en lumiere  introduit par des cristaux biréfringents, on trouve des 


polarisée, ’emploi de la lumiére blanche et Vobserva- 
tion des couleurs interférentielles permettent des 
mesures quantitatives. Si les deux ondes qui inter- 
ferent ont méme amplitude, les couleurs sont celles 
qwutilisait déja Newron pour la mesure des diffé- 
rences de marche [1]. Dans l’optique contemporaine, 
Vexemple le plus connu d’une telle mesure est la 
photoélasticimétrie, ou l’on obtient les courbes isochro- 
matiques correspondant aux retards de marche dus a 
la biréfringence du spécimen, mais la méthode est 
également importante pour les spécimens isotropes. 
Dans le dernier cas, le retard de marche optique est 


récents exemples de telles constructions dans le micros- 
cope interférentiel de Smiru [2] et de Nomarski [3], 
et dans les oculaires interférométriques de FrRan- 
con [4]. 

La méthode de FRANGON n/a pas, a notre connais- 
sance, été utilisée pour des mesures réelles dans le 
microscope, mais principalement pour rendre visibles 
des retards de marche entre différentes parties d’ob- 
jets transparents. L’un de nous [5] a utilisé cette 
méthode dans un montage macroscopique, pour me- 
surer la miscibilité et la diffusion mutuelle des liquides 
et étudier d’autres gradients du chemin optique. 


INGELSTAM 


E. 


Fic. 1. — Les couleurs de l’échelle de Newton, d’apres 
Baup-WricHtr [6], et quelques ellipses de discrimination 
dans la région sensible, d’aprés MAcAbAm [8], 10 fois agran- 
dies. Les retards de chemin optique sont donnés en microns 
sur la courbe. 


Dans Vouvrage cite, Ja lumiere spectralement 
pure du mercure a été utilisée, et non Ja lumiere 
blanche, essentiellement parce que la dispersion des 
détails pour la lumiére polarisée était génante et que, 
par consequent, la fidélité Wun équilibrage des cou- 
leurs n’était pas satisfaisante. 

En développant ces travaux et construisant un 
« diffusiometre » (1) extrémement simple, pour des 
liquides trés dilués, ces difficultés furent évitées. Nous 
nous sommes intéressés a la recherche de la précision 
avec laquelle on peut mesurer un retard du chemin 
optique comparant les couleurs de deux plages géo- 
métriquement semblables, l’une de ces plages formant 
Pimage du spécimen, l’autre constituant une plage de 
référence. L’intérét de cette recherche ne doit pas étre 
limité a cet appareil, les conclusions sont générales et 
la géométrie examinée ici est applicable 4 beaucoup 
@autres appareils interférométriques. 

Nous avons étudié particuliérement la teinte sen- 
sible du premier ordre et les teintes voisines de celle-ci, 
c’est-a-dire des couleurs rouges correspondant a une 
différence de chemin optique de 4500 A environ (on 
suppose l'appareil tel que I’échelle soit a centre noir), 
la teinte sensible a 5 700 A et les couleurs bleu-verda- 
tre 4 7 000 A environ. Baup et Wrrcurt [6] ont cal- 
culé les coordonnées trichromatiques et leurs résultats 
sont donnés sur la figure 1. En principe, on pourrait 
calculer les limites de discrimination effectives a 
partir des données fournies par Wricur [7] et par 
MacApam [8] qui couvrent tout le triangle des cou- 
leurs ; sur la figure 1, un certain nombre des ellipses de 
discrimination trouvées par MacApam sont tracées. 


(1) Pour la mesure de la diffusion d’un liquide dans un 
autre. 


ET 


[ Orr. Acta 
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On trouve que, pour le milieu de la teinte sensible a 
5 600-5 800 A environ, la fidélité d’une mise au point 
des teintes soigneusement faite doit étre de 10 A envi- 
ron. Pourtant, une étude directe nous a semblé dési- 
rable avec un appareil qui fait suivre aux couleurs 
échelle de Newton le long de la boucle de la figure 1. 
Une autre différence entre les mesures citées et les 
notres est que celles-ci étaient effectuées a luminance 
constante, alors que les teintes de ’échelle de NEwron 
changent constamment de luminance. I] semble néan- 
moins qwil soit raisonnable de limiter notre étude au 
domaine de la teinte sensible du premier ordre, parce 
qu’on voit directement d’apres le tracé que les circons- 
tances sont notablement moins bonnes pour les autres 
couleurs ; de méme, pour la teinte sensible du second 
ordre ou on ne peut prévoir qu’une précision d’un tiers 
de celle du premier ordre. 


Appareil. — L,/appareil utilisé, représente schema- 
tiquement sur la figure 2, est la modification dun 
dispositif antérieur [5]. La lampe a une température 
de couleur de 2 800° K. L’objectif O est formé de deux 
lentilles achromatiques trés soigneusement polies, qui 
font ordinairement les deux parois de la chambre de 
diffusion. La lentille achromatique L donne Vimage 
d’un plan voisin de O dans Je plan objet M d’un mi- 
croscope dans le champ duquel se trouvent les plages 
a mesurer. S est la plaque différentielle, constituée par 
deux plaques de quartz montées comme dans le pola- 
riscope de Savart. Elle est réglable, d’une part autour 
de Paxe optique, ce qui permet de régler la direction de 
biréfringence perpendiculairement a la fente d’en- 
trée E, et d’autre part, autour d’un axe horizontal 
paralléle a la plaque, ce qui permet de choisir avec 
precision une valeur quelconque de la différence de 
marche des rayons ordinaires et extraordinaires intro- 
duits par la plaque. Les deux polariseurs P, et P, 
sont toujours maintenus croisés et sous un angle de 
45° par rapport aux axes principaux des lames de 
quartz. Le quotient de l’intensité pour les polariseurs 
croisés et paralléles utilisés, est de 9.10-3 dans la 
région bleue et de 1,3.10-* dans le domaine rouge du 
spectre. 


_ Fic, 2, — Dispositif de mesure : E fente d’entrée ; O len- 
ree EDM L lentille formant limage d’un plan O en 
I, plan objet d’un microscope ; P, et P, polariseurs ; S pla 

diflérentieile, sa a ad is ‘al 
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Fie. 3. — Les coins de verre A et B suivant les lentilles O 
et les échelles circulaires. 


Dans le plan indiqué par la lettre O sur la figure 3, 
on place deux coins de verre A et B, qui doivent servir 
pour la mise au point des différentes couleurs. La me- 
sure est faite par leur rotation ; le coin A, qui doit étre 
reglé avec precision, est muni d’une échelle divisée en 
dixiemes de degrés, alors que B, qui doit étre orienté 
dans diverses directions, est muni d'une échelle plus 
simple. 

Si un coin, dont langle est e/y (voir fig. 3) et indice 
de réfraction est n, est tourné de l’angle y par rapport 
a la position de la figure, ot le gradient introduit par 
la plaque différentielle est maximum, le gradient du 
chemin optique dans le champ prend la valeur sui- 
vante 


ou b est le déplacement du a la biréfringence, et g le 
grossissement (ici inférieur a Punite, voir fig. 2). En 
vue de couvrir le domaine a étudier, quelques coins 
ont été fabriqués, leur pente e/y et indice de réfrac- 
tion n ont été mesureés. Pour le coin A (coin de mesure), 
le facteur du deuxieme terme de l’équation a une va- 
leur égale a 3.010 A. Ceci implique que, partant de la 
différence de phase réglée avec la plaque S (y = 90°), 
on ajoute ou on soustrait une fraction de ce montant 
en changeant l’angle ; pour des raisons de précision, 
on ne peut pas dépasser un écart de 70° de chaque 
cété. 

Pour éviter toute lumiere parasite notable, on 
masque entiérement l’ouverture a O, sauf deux carrés 
‘A et B (fig. 3). I] est nécessaire, pour la comparaison 
des couleurs des champs A et B qwils soient limités 
par une ligne droite ; a cet effet, l'image de B est jux- 
taposée a celle de A a l’aide d’un coin qui nest pas 
visible sur la figure. La perte de lumiere par réflexion 
est compensée par un verre plan de méme épaisseur 
interposé sur le trajet A. Le champ de l’oculaire dans 
le plan M peut étre limité en hauteur, et cette hauteur 
est choisie de facon qu’elle corresponde aux circons- 
tances des mesures. Les dimensions angulaires des 
plages A et B sont mises en évidence sur la figure 4. La 
luminance du champ est 3,0 nits dans les maxima pour 
la couleur jaune ; la lumiére parasite ne dépasse pas 
1/200 de cette valeur. 


SENSIBILITE DIFFERENTIELLE DE 


L’ GIL AUX TEINTES SENSIBLES 141 


Mesures et résultats. — On cherche d’abord, a l'aide 
de la teinte sensible, la position de la plaque A qui ne 
produit pas de gradient (y = 90°), et l’on introduit 
sur une partie des lentilles O le coin B qui a une pente 
plus grande que A. Par rotation de la plaque B on 
change la couleur du domaine B. Pour chacune de ces 
positions, on fait 10 mises au point successives de A, 
égalisant leurs couleurs. Comme mesure de la sensi- 
bilité différentielle aux couleurs, on prend la moyenne 
des écarts de ces positions par rapport a leur moyenne 
arithmeétique. Il est préférable de prendre comme 
mesure les écarts qui montrent le degré de reproduc- 
tibilté plutét que de tenter la mise au point d’une 
différence de couleur juste perceptible [8]. En géné- 
ral, on prend la valeur moyenne quadratique des 
écarts [8] conformément a la théorie statistique des 
erreurs. Cependant, quand il s’agit d’un petit nombre 
d’observations, il est tout aussi utile de prendre la 
moyenne simple. Nous avons de méme, pour une par- 
tie des observations, calculé la moyenne quadratique 
des écarts, et nous avons trouve que les résultats sont 
plus élevés @environ 35% étant donné les variations 
accidentelles. 

Les résultats de toutes les valeurs trouvées par deux 
observateurs sont portés sur la figure 5 (on a pris 
h, = 3°,2). L’écart moyen des angles est traduit en 
écart de chemin optique au moyen de l’équation 


} id AS) 3 O10 Rein | sdae sl = 


Pour fixer l’échelle des 4, on a cherché le premier 
minimum des interférences pour la raie 5461 A du 
mercure (obtenue a l'aide d’une lampe a mercure sui- 
vie d’un monochromateur). La plaque B était réglée 
de facon qu’on obtint une certaine intensité dans son 
champ interférentiel, c’est-a-dire que la différence de 
chemin soit un peu supérieure ou inférieure A 5 461 A. 
Ensuite, on a réglé la plaque A jusqu’a égalité photo- 
métrique avec le champ B dans les deux positions 
possibles : la moyenne de ces positions donne la valeur 
a l’échelle de A correspondant a la différence de che- 
min de 5461 A. Le réglage était fait avant et apres 
chaque série de mesures, ce qui permet de controler la 


microscope 
(la hauteur h est réglable). 


Fic. 4. — Le champ du 
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pour les teintes autour de celle sensible du premier ordre. 


stabilité de Vappareil (principalement celle de la 
plaque S qui est sensible aux changements de l’angle). 


Résultats. — Les mesures de la figure 5 ont été répé- 
tées plusieurs fois et, comme on le voit, il n’y a pas de 
variation dans la capacité de discrimination des deux 
observateurs. Le résultat le plus intéressant est la 
valeur minimum 4 la teinte sensible, ayant la valeur 
de 7 A environ. En outre, on trouve des valeurs consi- 
dérablement plus grandes pour les teintes voisines. 
Pour le rouge a 4500 A ov la courbe de cette couleur 
est paralléle a un coté du triangle de couleur, la valeur 
correspondante est 20 A, et pour le vert-bleudtre a 
6500 A, on a 25 A environ. La courbe est nettement 
dissymétrique, de sorte que Vincertitude est pius 
grande pour les teintes bleu-verdatre que pour les 
rouge-jaundtre, conformément aux conditions con- 
nues pour les couleurs pures et pour les zones corres- 
pondantes dans le diagramme trichromatique. En 
avancant de plus en plus au voisinage de la teinte sen- 
sible, on met au point aussi sur les intensités et non 
seulement sur les teintes de ces couleurs peu satu- 
rées [6]. Un petit nombre de mesures a été fait sur la 
teinte sensible du second ordre, l’écart moyen calculé 
de la méme fagon étant 30 A, c’est-a-dire trois fois 
plus grand que dans le premier ordre. 

Lorsqu’on mesure les coefficients de diffusion mu- 


tuelle des liquides, le champ n’est pas tout a fait uni- 
forme en phase, mais on mesure une région correspon- 
dant au maximum d’une courbe de Gauss [5], il est 
nécessaire de pouvoir limiter le champ de facon que la 
variation de phase soit négligeable. Pour cela, on a 
fait des études semblables avec une zone étroite 
(hk. = 0°,7) ce qui correspond a l’extension de la biré- 
fringence propre de la plaque S. 

Dans ces conditions, le coefficient de diffusion peut 
étre mesuré avec une bonne précision, en appliquant 
une petite correction [5]. En comparant les résultats 
obtenus avec la bande étroite, h,, a ceux que lon 
obtient avec h, pour toutes les valeurs entre 5 500 et 
5900 A du retard, on trouve une augmentation des 
valeurs. La moyenne des écarts pour h, étant pour un 
des observateurs, 7,1 A, a augmenté pour A, jusqu’a 
10,7 A, soit a peu pres 50%. A la teinte sensible, 
l’écart n’excéde certainement pas 12 A pour les zones 
étroites. 


Conclusions. — Les résultats nous autorisent a dire 
que l’erreur maximum obtenue dans une série de 5 a 
10 mesures reste de l’ordre du double de la somme des 
écarts moyens obtenus ici. Cela indique des préci- 
sions de 

20 A pour la teinte sensible, 

50 A pour les teintes voisines mentionnées, 
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Ces résultats doivent étre en bon accord avec ceux 
que l’on obtient par des mesures soigneuses sur des 
champs interférométriques de phase uniforme. Si 
toute l’échelle de Newron doit étre utilisée dans le 
domaine des trois premiéres périodes, on ne peut espé- 
rer qu'une précision de 100 a 200 A quand on estime 
le retard de chemin optique de l’égalité de deux cou- 
leurs, ou quand on essaie de suivre la limite entre deux 
couleurs voisines, si la phase du champ varie. A notre 
avis, il est absolument nécessaire @avoir un champ 
de référence pour égaliser une couleur; l’estimation 
« absolue » d’une teinte dépendrait trop des influences 
subjectives. 

Une suggestion trés importante est fournie par ces 
résultats. Si Pon veut effectuer des mesures précises 
avec ces types d’interférences, on doit essayer d’uti- 
liser la teinte sensible du premier ordre. Pour une autre 
difference de phase quelconque, cela peut se faire au 
moyen d’un compensateur convenablement construit. 
Le compensateur est réglé jusqu’a égalité avec un 
champ de références contenant la teinte sensible. 
Avec cette méthode, le compensateur constitue le 
dispositif de mesure et la mise au point soigneuse sur 
la teinte sensible constitue une méthode de zéro. 

En ce qui concerne l'appareil actuel, qui mesure les 
gradients directement, le compensateur est particu- 


SENSIBILITE DIFFERENTIELLE DE L’(IL AUX TEINTES SENSIBLES 


143 


liérement simple, parce qu’il a la forme de la plaque 
tournante A qui est placée devant le gradient dans le 
champ A a mesurer, dont on lit la position sur une 
échelle précise. Chaque gradient transformé en diffé- 
rence de phase a l’intérieur de la gamme du compen- 
sateur (par exemple jusqu’a 20 000 A), est donc acces- 
sible pour la mesure avec la méme limite d’erreur 
(20 A). On corrige facilement la différence de disper- 
sion des deux milieux. Cela veut dire que des coeffi- 
cients de diffusion des solutions extrémement diluées, 
sont mesurés avec une précision aussi grande, c’est-a- 
dire a quelques milliemes pres, que celle des méthodes 
dinterférences ordinaires utilisant une différence de 
marche de plusieurs longueurs d’onde fournie par les 
solutions plus concentrées. 
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INFORMATION 


Colloque International de spectroscopie 
Amsterdam, 14-19 mai 1956 


Le sixiéme Colloquium Spestroscopicum Interna- 
tionale sera organisé 4 Amsterdam pendant la semaine 
qui précéde la Pentecdte (14-19 mai 1956) par les soins 
des Sociétés Royales Néerlandaises de Physique et de 
Chimie. 

De plus amples détails concernant ce colloque, ainsi 
que les formulaires d’inscription ne seront envoyes 
qu’a ceux qui, en communiquant leur nom et leur 
adresse au Bureau du Colloque, auront manifesté leur 
intérét a cette réunion. 

Le Colloque traitera de la spectroscopie d’émission, 
de la spectroscopie d’absorption, principalement en 
vue de l’application a l’analyse. 


Les personnes qui désireraient présenter une commu- 
nication doivent se faire connaitre le plus tét possible. 

Pour tous renseignements, s’adresser au 

Prof. Dr F. Fresse, Secrétaire de la Commission du 
Colloque International de Spectroscopie. 

Laboratorium voor Analytische Chemie, Nieuwe 
Achtergracht 25, Amsterdam C4. 


Plan de la réunion. — 1. Spectroscopie d’émission : 
a) la photométrie de la flamme. 6) l’arc électrique. 
c) l’étincelle. d) la quantométrie. e) linstrumentation. 
— II. Spectroscopie d’absorption : /) la spectroscopie 
infrarouge et la spectroscopie RAMAN. g) la spectro- 
scopie ultraviolette et la spectroscopie de fluorescence. 
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Application of the binocular matching method to the study 
of the subjective appearance of surface colours 


Elaine G. T. WaASSEF 
Institute of Education for Girls, Ibrahim University, Cairo, Egypt. * 


SumMARY. — A binocular matching technique has been developed for application in studies of the effects of visual adaptation on the 


appearance of surface colours. 


The assumptions underlying the technique were considered. A factorial experiment was designed to examine the relative impor- 
lance of the influences of changing the illuminant and the background on the hue of some red, green and purple samples of diffe- 
rent levels of value and chroma, and to assess the interactions of these factors, Statistical analysis of the data showed that the 
standard error of a single match was about two hue steps ; and that it was possible, without repeating any matches, to evaluate 
the influences of the above mentioned factors within close tolerances. 


Sommaire, — Une méthode de comparaison binoculaire a été mise au point pour Vapplication aux études de Veffet de Vadaptation 


visuelle sur Vapparence des couleurs de surface. 


Les hypotheses sur lesquelles est basée la technique expérimentale ont élé éludiées. Une expérience a été concue pour déterminer 
les importances relatives des changements d’éclairage ou de fond sur la teinte de quelques échantillons rouge, vert et pourpre de 
différentes saturations ou couleurs, ainsi que pour délerminer les interactions entre ces facteurs. L’ analyse statistique des résul- 
tats a montré que Verreur moyenne dun équilibrage élait de Vordre de deux échelons d2 couleur et qzwil est possible, sans répéter 
aucun équilibrage d@’évaluer les influences des facteurs ci-dessus avec une bonne précision. 


ZUSAMMENFASSUNG. — Hine binokulare Vergleichsmethode wird entwickelt und auf die Untersuchung angewendet, welchen Einjluss 
die visuelle Adaptation auj die Erkennbarkeit von Oberfldchenfarben hat. Die Voraussetzungen, auf denen sich diese Methote 
aufbaut, werden gepriift. Ein besonderer Versuch zeigt die relative Bedeutung von Anderungen der Helligk*it und des Hinter- 
grundes auf die Schattierung einiger roter, griiner und purpurfarbenen Prob n verschiedener Wertstufen und Farben. Man kann 
die Wechselwirkungen dieser Ein/liisse abschdtzen. Eine statistische Analyse der Erg?bnisse zeigt, dass der mittlere Fehler eines 
einzelnen Vergleichs ungefdhr zwei Farbstufen betrdgt. Es ist méglich, ohne Wiederholung der Vergleiche die Finfltisse der 


erwdhnten Faktoren innerhalb enger Grenzen zu ermitteln. 


1. Introduction. — The efficiency of the binocular- 
matching method as a tool of investigation in vision 
was established by Wricurt [1]in 1934, but its appli- 
cation to the study of the appearance of surface co- 
lours seems to have been first attempted by Bouma 
and KruitHorr [2] in 1943. The method has the 
advantage of yielding more consistent observations 
than the memory-matching methods which have 
been widely used in this field [3, 4, 5]. That experi- 
menters have not made much use of the binocular 
method may be due to the feeling that memory mat- 
ching is the natural way of recognizing colour changes 
in everyday life. It is also possible that opinion was 
influenced by Bouma & KrurrHorr’s [2] own con- 
clusion that the binocular method could not replace 
the memory method in view of the discrepancies bet- 
ween the results which they obtained by the two 
methods. In the writer’s opinion the discrepancies 
recorded by Bouma and Krurrnorr are of doubtful 
significance since the comparison was made between 
the data after the memory observations had been 
smoothed by means of the moving-average formula 
which, according to KrnpaLt [7], may generate 
trends of no physical origin. 

The work described in this paper was an attempt 
to re-assess the power of the binocular method in 
connection with studies of the subjective appearance 


* Formerly, Technical Optics Section, Imperial 


College, 
London. 


of surface colours. For this purpose a binocular viewed 
was constructed, and the experiment was planned 
and analysed according to modern statistical views. 


2. Preliminary tests. — The following observations 
were made on Wriacut’s colorimeter [8] to test (a) the 
possible effect of the state.of adaptation of one eye 
on colour appraisal by the other eye, and (b) whether 
the observer’s eyes differed in their colour perception, 


Test (a) : The relative luminous-efficiency values 
of the right eye were measured for a number of wave- 
lengths when the left eye was dark adapted in one 
case and light adapted in the other. The values were 
smaller in the latter case, but the differences did not 
reach the 5 per cent significance level except at a few 
wavelengths. Nevertheless, the combination-of-pro- 
babilities test substantiated the indication of the per- 
sistent sign that the state of adaptation of one eye 
had an effect on the other eye. While this agreed 
with the conclusion of AguiLar and Soxis [9], it 
showed that the effect was not large enough to in- 
fluence the results drawn from observations with the 
binocular viewer. 


Test (b) : Each eye saw a separate field illuminated 
by a monochromatic radiation drawn from an inde- 
pendent spectrum. The recorded differences between 
the wavelengths of the illuminating radiations when 
the two fields were matched indicated the presence of 
small differences in the colour perceptions of the two 
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eyes. The precaution was therefore taken to match 
each reference sample once with the left eye and once 
with the right eye. 


3. Apparatus and observational procedure. — The 
binocular viewer (fig. 1) consisted of two adjacent 


Fic. 1. — Binocular Viewer (not to scale) 


compartments approximately 36” x 36" x 36” each, 
separated by a light-tight partition. The inside of the 
two compartments was painted black. Each com- 
partment was illuminated by two pairs of lamps fixed 
in directions making angles of 45° at the centre of 
the side wall where the sample was to be placed. The 
lamps gave a fairly uniform illumination over the 
size of the field used in the experiment. The light 
reflected from each sample and background entered a 
right-angle prism where it was reflected through 90° 
and passed through two openings in the front face of 
the viewing box. The distance between the openings 
was adjustable to suit the observer’s interocular dis- 
tance. The symmetrical design of the compartments 
made it easy to interchange the samples and condi- 
tions of viewing. 

The lamps used in the experiment were gas-filled 
incandescent pearl lamps matched to a standard lamp 
calibrated at the National Physical Laboratory to give 
the colour temperature of the CIE Standard [lumi- 
nant A (about 2 850° K). Illuminants B and C were 
obtained from combinations of the standardized 
lamps and ‘“‘ Chance OBY ” glass filters. The colour 
temperature of the lamps was checked after using 
them for thirty hours, but very slight changes were 
found. 

Colour samples from the 1950 edition of the Mun- 
‘seLL Atlas were used throughout the investigation. It 
contained forty charts of constant hue, at intervals 
of 2 1/2 hue steps (The hue circle is divided into 
100 steps starting from and ending at red). Each chart 
contained samples of the same hue at intervals of 
one step in value (1, 2,..., 9) and two steps in chroma 
(0, 2, 4...). Some samples of chroma 1 and 3 were also 
available in the ten major hues. 

The observer sat in a dark room and adjusted the 
distance between the openings in the front wall of the 
viewer to suit his interocular distance. Each eye then 
viewed a field subtending an angle of 22° x 17° with 
an opening of 1° x 1° 20’ in the middle of the field. 
A final adjustment of the separation between the two 
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fields was made by small angular rotations of the 
prisms to avoid fusing the fields. 

Each session of observations was preceded by a 
period of five minutes light adaptation of the eyes to 
the conditions in the respective compartments. The 
reference sample was then pressed against the wall 
of one of the compartments on the small opening to 
reduce the shadow due to the space between the back- 
ground and the sample. To make the match the ope- 
rator presented to the observer in the other compart- 
ment a MuNsEeLL sample from the Atlas charts of 
constant hue. The observer was asked to compare it 
with the reference sample viewed by the other eye, 
and if it did not match he would ask for a lower or 
higher value or chroma or a specific change of hue. It 
was found that an attempt to find an approximate 
match of value and chroma from the first chart 
usually facilitated judgment of the required change 
of hue. If the required match was not found among 
the samples, it was estimated by interpolation in all 
three attributes. In some cases, however, extrapola- 
tion was inevitable. After the observer had decided 
on the match, the operator recorded it in terms of the 
Munsell notation. Another reference sample was then 
presented to the observer and the process of mat- 
ching repeated. About ten reference samples were 
matched in each sitting. The samples were presented 
in a random succession, and it was also observed that 
no match was repeated in the same sitting. Sometimes 
two sittings were made in one day. 

It is important to note that the observer could not 
fixate on both fields at the same time, i. e. the method 
was not of simultaneous comparison. It was in fact a 
very quick memory matching where the standard 
was always available for comparison. To avoid adap- 
tation effects due to the samples themselves the obser- 
ver did not look at either for any length of time, and 
the samples were not left in the openings during the 
five-minute adaptation periods. 


4. Design of the experiment. — The appearance 
of a surface colour depends essentially on two factors : 
the light reflected from the sample, which is physi- 
cally measurable, and the state of adaptation of the 
eye. Earlier work on adaptation using aperture co- 
lours [8, 11, 12] and surface colours [2, 3, 4, 5, 10] 
showed that the state of adaptation of the eye depends 
on several factors : the type and intensity of the adap- 
ting radiation, reflectance and type of objects in the 
field of view (i. e. the background), and the size and 
type of the adapting field. The duration of the expo- 
sure of the eye to the adapting radiation and test field 
also affects its state of adaptation. The response to 
these factors may differ from one person to another. 
Furthermore, their effects are not necessarily inde- 
pendent. For example, the influence of changing the 
illuminant on the appearance of the sample depends 
on the reflectance of the background. Such a pheno- 
menon is referred to as an interaction between fac- 
tors, in this case between the effect of the illuminant 
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TABLE 1 


Matches with white background 


[Orr. AcTA 


Left eye matches 


Right eye matches 


A C B A C B 

SOR 3/2 4 YR 3.8/2.5] 9 R Soyer (tay NAR 3.5/3 5) WAR 4/2 6R 3.8/2 oRP 4/1.5 
3/4 10 R 5.5/4 6.5 R 3/6.5 |8.5 R Spy fes | MOR 3.8/4 oR 3.5/6 6R 1.5/5 
3/6 19 R 3.9/8 6R 4/8 6R 1/8 OR 4/5 dR 4.5/7 2.0 R 4/6 
5/2 3.5 YR 5.5/2 7.0 R 6/3 2.9 YR 9/3 10R Deal Tos) IR 6/3 6.5 R 6/2 
5/4 Ia) WARE IZ 7.0 R Deo | Toa) IR 5.0/7 SO) LR 5.5/4.5} 2.5 R 5.5/6 6R 5.5/6 
5/6 10R 5/8 5.0 R 6/11 6R 5.5/10) | 720 5.2/7 2.0 R 5/10 1.5R 5.7/10 
7/2 10R 8/2 Ou CPIBS Is 7.8/2.0)| OX 8/3 10R 8/3.5 TR 8/5 
7/4 6R 7.0/7 6.5 R 8/7.5 |6R DP a LOREX 7.0/9 10 RP 7.2/8 3.0 R 7.5/8 
7/6 7.9 R 7/7 6R Heo) Oe Onix 7.5/10 |10 R 7.2/6 Zed EX 7.2/8.5| 4R 8/8 

5.0G 5/2 2.0 G 5.0/3 3.0 G 6.5/2 |4G 5.2/3 6G 5.5/2.5) 7.5 G 6/3 ONG: 5.8/2 
5/4 1.5G 5.5/6 Mea; DOO foro) Cx 5.8/6 3.0 G 5.8/5 |10G 5.5/4 |10G 5.5/4 
5/6 SOG 5.8/8 748) (G; 5/7 6G 6/8 2.0 G 6.2/5 8G 5.5/6 9G 6/4 
7/2 2.0 G& 8.5/5 5G 7.8/2 |5 G 7.8/3 2.5 G 8.5/4 2.9 BG 7.8/2.5) 7.5 G 7.8/2 
7/4 3.5 G 8/7 6G 7.8/6 |5G 7.8/6 4G 8/5 7.5 G 7.5/5 6.5 G 8/4 
7/6 2.0°G: 7.5/8 6G Hesy7h \\Mesy KG; 7/8 oG 8/6 1:0, G 7.2/6 |10G 7/5 

A ,(0) 12 3/2 one 3.5/1 6 RP ayceypile [ay SIR 3/2 2.5 P 4/1.5 ye 3.8/2 8.5 P 3.5/2 
3/4 We Ne 3.0/3.5] 6.5 RP 3/4 6 RP 3.2/3 2.5 P 4/2 LORE. 3.0/4 |10 P 3.5/3 
3/6 (3) 1 3.5/6 330 RP 3.0/6" 4 3/4.5 Dae 3.5/4 2.0 RP 3.5/4 ge 4/4 
5/2 ilioy 1? 5.0/1 ey UR 6/1.5 |5 RP Sy PHietay))aligay Vee 5.5/2 hod 6/2 TRONS 6/1 
5/4 De 5.5/4 4 RP 5/4 6 RE 5/3 Se) 1 5.0/3.5/10 P 3.0/0-0)| 9 P 6/3 
5/6 ale 9.0/6 _|10 RP 5/5 o RP 5/5 4P 5.2/6 |10 P 5/5 IER 5/4.5 
7/2 i) Ne? Seo) leo ones 8/2 6 RP 7.9/2 1.5 P 8/3 8.5 P 8/2 6 P 8/2 
7/4 by de? 8/5 3.0 RP 8/5 DRE 7.2/5 6P 8/4 10 P 7.5/6 IPR 7.9/6 
7/6 Tea) V2! 6.8/6 6 RP 7/9 2.5 RP WS ge 6.8/7 Ss) AP 5.5/6 Oe 1.0/0 

TaBLe 1 (contd.) 
Matches with black background 
Left eye matches Right eye matches 
A C B A G R 

OOM SIZ 4 RP 4.5/5 6R 4/4 oR 4/4 8.5 RP 4/4 6R 3.8/4 6R 3.8/4 
3/4 1.5 R ANG || a 4.5/8 5R 4/7 25R 4.2/6 4R 4/9 4R 4.5/6.5 
3/6 4R 42/10; 5.5 R 4/10 4R 4.2/9 15 Re 3.0/8 4R 4/9 4R 4/9 
5/2 Wd) 1 6.5/5 8.5 R 6/5 SeONky 6/6 LORE 5.8/4 6.5 R 6/5.5 4R 6/5 
5/4 PH Iae  Saay lle) || (fy 1a} 5.5/9 6R 5.0/7 7.09 RP 5.5/6 6.5 R 5/7 3R 6/10 
5/6 7.9 RP 5/12 6R D2 LOM mone 5.5/10 | 7.5 RP 5.5/10 | 6R 5/10 6.5 R 5/8 
7/2 ToD AP 8/6 10R 8/8 2.0 R 8/10 Oe 7.2/8 7.0 R 7T.2f 1 (pda 7.5/5 
7/A De OMEN eae | lc 7.5 R 7.2/10 | 6 R 72/10) )/10 P 7/8 6R 7/@ 3R 7/10 
7/6 ORE it 8.5 R 7.2/10| 5 R 8/10 2.5 RP 6.5/8 6.5 R 7/10 2.09 R 7/10 

HOG, ay 3.5B 6/3 1.5G 8/1.5 8.5 G 6.2/2 |10 BG 5.5/2 UGE 5.5/2.5] 6G /2 
5/4 4 BG 5.8/6 6G 5.8/3.5] 7.5 G 5.5/4 1.5 BG 5 5/4 3G aia! 5G 6/3. 
5/6 4B 5.5/7 6G 5.8/6 |10G 5.8/6 6G 9.9/5.5] 5G 5/7 5.5 G 6/5.5 
7/2 10 B_ 8/6 2.5 G 9/3 4G 9/3 2.5 B 7.5/4 |10 GY 8/2 10 GY 8/3 
7/4 3.0 B 7.2/8 5G 8.5/6 5G 7.5/6 8.5 B 7/4.5 4G 7.8/4.5) 5G 8/5 
7/6 7.9B 7/8 4G 9/8 ONG, 7.5/8 7.5BG 6.5/6 5G 7.5/5 5.5 G 8/6 

SAO) Ve YP) one 5/4 6.5R 4/2 10 RP 3.8/2 8.5 P 3.5/2 2.5 R 4/2 EN. 
3/4 4P 4.8/9 |10 RP 4/5 5 RP 4/5 ue 35/6 8.5 RP 3 3/6 9.5 RP Fi 
3/6 one 4.8/9 8.5 RP 4/6 4 RP 3.8/6 4P 4/9 5 RP 3.8/7 GIGS) hfe fas 
5/2 1.5 1 6.2/6 7.0 R 7/4 Laue! 7/6 5» P 5.8/4 7.5R 5.5/3 3.5 RP 6/3 
9/4 3.5 P 6.2/9 10 RP 6.2/5 5 RP 6/5 4P 5.5/5 3.5 R 5.8/6 2.5 RP 6/5 
5/6 5P 6/10 1 R 6/7 4 RP 6/6 aye 5.5/7 6.5 RP 6/8 2.5 RP 5.2/6 
7/2 4P 8/6 6.5 R 8/5 6.5 R 8/6 2ne 7/6 2.5 R 7.2/7 PY VED FSIS) 
7/4 2.9 P 7/10 3.5 R 7.2/7 eRe 8/7 2.5 P 4/6 4R 7/6 1.5 RE 8/8 
7/6 5P 7/9 10 RP SHAG) IRS 7 fsyAl(Oy || 4b 18) 7/8 ORE 7/9 IYI RNS 7s 


———————————ooo — 
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and the effect of background. As the possibility of 
such mutual dependence cannot be ruled out it be- 
comes imperative to study the possible interactions as 
well as the general effects of the factors. Since a detai- 
led study of each factor in this way would require a 
large amount of experimentation, it was thought ne- 
cessary in the first place to sort out the more impor- 
tant factors and interactions by means of an expe- 
riment planned according to the factorial scheme 
which has proved to be very efficient in similar types 
of problems in other fields of research. 

The factors covered by the experiment were 
(1) the type of illuminant, by experimenting on illu- 
minants A, B and C, at approximately 10 lumens per 
square foot, and (2) the reflectance of the background 
by using a white and a black background. Similar 
backgrounds were used for the reference and mat- 
ching samples. The reference samples were always 
illuminated by illuminant A, and it was found prac- 
tically more convenient to maintain its level of illu- 
mination at approximately 36 lumens per square foot. 

In view of the slight difference in the colour per- 
ception of the observer’s eyes (Section 2) all the 
matches were repeated after interchanging the vie- 
wing conditions in the two compartments of the vie- 
wer. The difference between the eyes was treated as a 
separate factor when analysing the data. 

The three major hues 5.0 R, 5.0 G and 5.0 P of the 
MunseE Lt Colour Atlas were experimented upon. Nine 
samples from each of the red and the purple hues, and 
six samples from the green were observed, correspon- 
ding to the nine different combinations of the chro- 
mas /2, /4, /6 and the values 3/, 5/ and 7/ with the 
exception of values 3/ in the green because the sample 
5.0 G 3/6 was not available. This choice made it pos- 
sible to investigate the influence of the various illu- 
minant-background combinations on, say, the hue at 
various levels of value and chroma. 

288 matches were made, but no match was repeated 
under the same conditions of viewing except for occa- 
sional checks on the consistency of the observations. 


5. Statistical significance of the factors. — The 
matches are listed in Table 1. The analysis was confined 
to the changes of hue. The hue shift, which is the hue 
. of the matching sample minus the hue of the standard 
sample, is tabulated in Table 2 in terms of MUNSELL 
steps of hue. 

A few general trends are discernible in Table 2, but 
specific results can be obtained only by statistical ana- 
lysis of the data. The method of the analysis of va- 
riance was therefore applied to investigate the varia- 
tion pertaining to each hue and background. The 
purpose of the analysis was to determine the relative 
importance of the factors : illuminant, chroma, value 
and the eye-effect ; and at the same time to assess 
their interactions and to estimate with precision the 
error of a single match which is a criterion of the repro- 
ducibility of results obtained by binocular match- 
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ing. In the light of this analysis we shall evaluate the 
confidence limits of the matching hues. 

The hypothesis of colour constancy served well as 
the null-hypothesis of the analysis of variance. If it 


TABLE 2 (a) 


Difference, in MunseLt hue steps, between the hue of the 
matching sample, seen by the left eye against the white 
bacground, and the hue of the standard sample seen by 
the right eye against an identical background. 


IIluminant A Iluminant © Iluminant B 

R G P R G Pp R G P 
3/2 9.0 0 4.0 11.0! 10.0 30.0 
3/4 5.0 Zoi a0) SII) ae) 11.0 
3/6 9.0 HEH (0) 8.5 1.0 19.0 

| 

0/2 | 8.5 |—2.5/—3.5| 2.6 |—1.5) 10.0] 7.5/—1.0] 10.0 
5/4 (en) |i) 0 Pro) Peta) 9.0 2.5;—1.5] 11.0 
5/6 5.0 |—1.5 0 0.5 ZrO eon) 1.0 HCO) WOM) 
WE | sell) Dees] 0) 1.0 0 4.0 || 0) 11.0 
7/4 Ce eee 0) ILE) 1.0 8.5 0) 0) 10.0 
7/6 2.0: |—=2.0 Zao bal) WOK Tal) 5) Te 


TABLE 2 (b) 


Difference, in MuNSELL hue steps, between the hue of the 
matching sample, seen by the right eye against the while 
background, and the hue of the standard sample seen by 
the left eye against an identical background. 


IHluminant A Iluminant C [luminant B 

R G P R G P R G P 
3/2 10.0 —2.5 1.0 0) —6.0 OrO) 
3/4 0 —2.5| 0 nd) 1.0 5.0 
3/6 1.0 —3.0) 0 7.5 |—2.5 4.0 
Ay) 5.0} 1.0/—3.5] 2.5] 2.5 | 2.5 alia 240) | Piss 
5/4 4,0|—1.5|—1.5|—2.5| 5.0 | 5.0 NAO tO | AAO) 
5/6 PS =P 5) 1) = 3) SO) || EO) seas) 20) ) (ow) 
7/2 4.0|—2.5|—3.5| 5.0) 7.5 | 3.5 |—4.0) 2.5 | 1.0 
7/4 5.0/—1.0) 1.0/—5.0} 2.5 | 5.0 |—1.5)] 1.5 | 6.0 
7/6 5.0} O 3.0] —2 2.5 | 4.0 |—1.0) 5.0 | 4.0 

i) 


TABLE 2 (c) 


Difference, in MUNSELL hue steps, between the hue of the 
matching sample, seen by the leff eye against the black 
background, and the hue of the standard sample seen by 
the right eye against an identical background. 


IWuminant A Iluminant C Iluminant B 

R G Pp R G Pp R G eg 
3/2 |—11.0 0) Pay (0) 15.0 
3/4 |— 3.5 —1.0 1-0) ee 10.0 
3/6 |— 1.0 0 ).5 13.5)—1.0 9.0 
5/2 |—17.5| 18.5|—3.5| 3.5|—3.5| 22.5/—1.5| 3.5] 11.0 
5/4 |—11.5| 9.0|/—1.5 il 35) LOH) abs 0) 1.0} 2.5) 10.0 
5/6 |— 7.5 9.0} O LO} 150) 1620) 50 5.0; 9.0 
7/2 |\—17.5| 25.0/—1.0} 5.0)—2.5) 21.5|—2.5|—1.0| 21.5 
7/4 |\—12.5| 18.5|—2.5| 2.5) 0 18.5 TX) (0) 14.0 
7/6 |—13.5] 22.5) 0 13.5)/—1.0] 15.0} 0 AAD "7385 
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TABLE 2 (d) 

Difference, in MUNSELL hue steps, between the hue of the mat- 
ching sample, seen by the right eye against the black back- 
eround, and the hue of the standard sample seen by the 
left eye against an identical background. 


[luminant A Tluminant C [lluminant 

R G 1p R G P R CG Pp. 
3/2 |— 6.5 oD) 1.0 INS) 1.0; 6.0 
3/4 |— 2.5 O |—1.0 13.5|/—1.0 TES) 
3/6 |— 3.5 —1.0/—1.0 i10.0}—1.0 6.5 
5/2 |—15.0] 15.0} 0 1.5|—3.0) 22.5)—1.0 TRONS SS 
5/4 |— 7.5 6.5|—1.0 1.5|—2.0|} 18.5|—2.0| 0O Ws 
5/6 |— 7.5 lea) |e) Ol INATENSs Toss) (esi 7a 
7/2 |\—16.0| 17.5|—3.0] 2.5|—5.0} 17.5 1.0|—5.0| 7.5 
7/4 |\—15.0|} 13.5|—2.5 1.0)—1.0} 19.0)/—2.0| 0 Oso 
7/6 |—12.5| 12.5)}—1.0 1.0} O 15.0/—2.5| 0.5! 6.5 


were true all the hue shifts would be nil within the 
limits of experimental error. The analytical procedure 
of the analysis of variance is to split up the sum of the 
squares of the deviations of each shift from the mean 
of all the shifts into (a) components corresponding to 
the factors which cause (or are suspected to cause) 
variation, (b) other components corresponding to the 
interactions of these factors and (c) a residual from 
which the experimental error can be estimated. Each 
sum of squares is divided by the respective number 
of degrees of freedom, and the quotient is the ‘‘ mean 
square ”. When the mean square corresponding to an 
interaction proves to be significantly larger than 
the mean square of the experimental error the infe- 
rence is that the influences of these factors are not 
independent, i. e. the effect of one factor depends on 
the level of the other factor. When a factor has a 
general effect, i. e. when its influence appears at all 
the levels of the other factors, its mean square is 
expected to be significantly larger than the mean 
square corresponding to its interaction with them. The 
standard experimental error is the square root of the 
mean square of the residual and has its number of 
degrees of freedom. 

The analysis of variance of the data was carried out 
along the lines sketched above, the computing scheme 
being as described by Brown.ee [13]. The details of 
the calculations are omitted, and an example of the 
analysis is given in Table 3 where the letters I, C, V, 
and E stand respectively for the influences of illumi- 
nant, chroma, value and the difference between the 
estimations of the shift by the two eyes. FisHer’s 
z-test was applied to each mean square, and the resul- 
ting levels of significance are indicated in Table 4 by 
one, two or three crosses depending on whether the 
level of significance is 5 per cent, 1 per cent or 0.1 per 
cent, respectively. 

The persistent significance of the EI interaction 
indicates that the effect of the difference between the 
colour perceptions of the eyes changed from one illu- 
minant to the other. Since neither EV nor EC was 
significant one infers that the eye-effect did not de- 
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pend on the value or the chroma of the observed 
sample — within the ranges covered by the ex- 
periment (?). 

The significance of the interactions IC and IV in the 
larger proportion of the investigated viewing condi- 
tions shows that the hue shifts associated with the 
change of illuminant (A, B or C) varied with the chroma 
and the value of the sample. The fluctuations of the 


TABLE 3 
Analysis of Variance of the Hue Shifts of the Red Samples 


Source Sum of Squares Mean Squares 
of dat 
Variation Black White Black White 
ss 1 46 92.04 -46 92.04 
I 2 1 384.62] 221.36 692.31 110.68 
V 2 64.12 18.69 32.06 9.34 
C 2 27.37 95.44 13.68 47.72 
[Dix 2 12.40 53.08 6.20 26.54 
ESaVi 2 aay 4.86 1.90 2.43 
EXC 2 9.99 appa la | 4.80 6.06 
DV 4 253.66 12.03 63.42 3.01 
Isa 4 116.07 13.36 29.02 3.34 
Nea, 4 10.32 ale ital 2.58 2.78 
Residual 28 47.08) 147.29 1.68 5.26 
Total 53 1 929.48| 681.37 
TABLE 4 
Significance Levels of the Effects 
Source Red Samples Green Samples Purple Samples 
of 
Variation | Black | White Black | White | Black | White 
E a — we 
I 2 — ae 
Vv = — = 
C — = = 
+++ = — 
Exl + =e ee AE 
++ ~ ++ 
EX Vi — — = 
- — + 
EXC = = mt 
rxv 9 |e +++ = 
-- = i 
PG ee +++ ++ 
Vee — — = 
Mean ae 
Square 7) 4.09 5.35 
Error 4.56 1.92 9.60 


(1) In fact EV was significant on the 5 


— possibly a chance effect, 


per cent level once 
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shifts with the factors V and C, as well as with K, 
were in fact so large that the change of illuminant 
did not demonstrate a significant general effect 
except for purple and possibly red samples with a 
black background. The dependence of the effect of 
changing the illuminant with the chroma of the 


sample was more pronounced than its variation with 
the value. 


6. Numerical evaluation of the effects. — In the 
hight of the above discussion of the statistical signifi- 
cance of the factors and their interactions the ranges 
(of hue) within which the matches should be expected 
to he can be calculated. 

In each case three quantities were calculated : the 
mean value of the shift, the standard error of this 
mean and its 95 per cent confidence limits. In general 
the average of the matches observed by the two eyes 
was taken as the match which would have been found 
under the normal conditions of binocular viewing. 
Also, in the cases where the hue shifts did not vary 
significantly with the level of a factor the mean of 
the shifts recorded at the various levels of this factor 
gave the most probable estimate of the shift. Advan- 
tage was thus taken of the presence of insignificant 
effects to increase the precision of assessing the more 
influential factors. 

Consider for example the red samples against a 
white background. In this case V and its interactions 
were not significant. The mean of the shifts recorded 
at the three levels of the value were therefore used. 
Taking, in addition, the mean of the matches observed 
by the two eyes the relative data could be condensed 
and are given in the following table where each item 
is the mean of six measurements of the hue shift. 


A B C 
Ips 6.9 DP BP 
[4 4.4 iL = 
[6 4.0 — .7 — 4 


The standard deviation of a single observation is 
2.14 hue step (from Table 4). Hence the standard 
in the table is 2.14/(6)2 = .87 
hue step. The value of Student’s t-statistic for 
5 degrees of freedom (number of observations mi- 
nus 1) at the 95 per cent probability level is 2.0. 
. Hence the 95 per cent confidence limits of each item 
is equal to the mean + .87 2.0. It follows that when 
the 5.0 R samples of chroma 2 are matched under 
illuminant C we should expect that 95 per cent of 
the matches will have their hues within the range 
5.0 R + 3.2 + 1.74 = 8.2R + 1.74. 

The matches made with the two backgrounds were 
treated in this way with the following results : 

(a) Case of white background. — The confidence 
limits of the hues of the matching samples when mat- 
ching is carried out against a white background are 
given in Table 5. 

An examination of this table shows that : 

(1) There were no significant differences between 


deviation of each item 
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TABLE 5 


95 % Confidence limits of the hues of the matching samples 
in the case of the white background ‘based on the means 
of sets of six observations in the cases of the red and the 
purple samples, and on sets of four observations in the case 
of the green samples). 


A B ( 
50K (2 MONR 2 2079 R Aaa ph aan 
A 19.4 R 2 236 R 2 Sol4a Re Bee 
6 SOR 2 2045p = 954 6 eer 
BOG 2 (346-2 292646 = 22\7 Gees 
if \81-G 2 23160 6 25-207 Gee 
16 \S4G 2 PolsrG 40 99170 G aeeeeee 
50P 3 49P 2.32/71 RP + 32/00 RP eee 
5} |84 P= 4° 3.2/2.2 BP = 3108 Rpeeeree 
4] ISS P 4 896 RP + 3,211.0 RPeeeree 


the appearance of colours under illuminant B and 
their appearance under illuminant C. 

(2) The influence on hue of changing the level of 
illumination (which was exemplified in this experi- 
ment by matching under illuminant A of lower level 
of illumination) differed considerably from one hue to 
the other. The hue shifts were large in the case of red 
samples, smaller for green samples and insignificant 
in the case of purple samples. 

(3) The appearance of colours under the illumi- 
nants B and C differed significantly from their appea- 
rance under A; red and green samples looked more 
yellowish under B or C, while purple samples looked 
bluish. 

(b) Case of black background. — In view of the large 
number of significant factors and interactions in the 
cases of the red and the green samples the original 
data could be condensed only in as far as taking the 
means of observations made by both eyes. The means 
of the matching hues are given in Table 6 for all 
three colours. The 95 per cent confidence limits, based 
on six observations, are given in the case of the purple 
samples where V and its interactions were not signifi- 
cant. 

Comparison of Table 5 and Table 6 shows that : 

(a) The changes in the appearance of colours are 
in most cases greater when viewed against a black 
background than when the background is white. This 
may be explained by the greater influence of adapta- 
tion in the latter case due to the larger amount of 
light reflected by the background which reduces the 
effects of the colour rendering properties of the illu- 
minant. 

(b) Only in the case of black background have 
significant differences appeared between the hue 
shifts observed under illuminants B and C the spec- 
tral distributions of which are very close. 

(c) In some of the cases the change of background 
reversed the direction of hue shifts. 


7. Conclusions. — The preceding experiment de- 
monstrates the efficiency of binocular matching for 
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TABLE 6 
Mean of the matching hues in the case of the black background 


, : : : i : B C 
3/ 5/ 7/ 3/ 5/ 7/ 3/ 5/ 7/ 
5.0 R {2 6.2 RP reyes LE? 8.2 P 5.0 R 3.8 R 4.2R 6.0 R 7.5 R 8.8 R 
/4 2.0R Dron bale ies Leas 4.5 R 4.5 R 4.5R 5.0 R 6.5 R 6.8 R 
/6 2.8 R Wd IRE P20) Re 4.0 R 5.8 R 3.8 R 4.8 R 6.0 R 7.5 R 
AWE, 2 — 1.8B 6.2 B —= 7.2G 2.0G — 9.2 GY 1.2G 
/4 — 2.8 BG 1.0B = 6.2 G 9:0 1G — 4.5G 4.5G 
/6 a= 10.0 G 2.9 B — 7.8G 7.2 G — 5.5 G 4.5G 
5.0 P 2 43P + 2.4 6:6 RP + 2.4 oo Rae 2:4 
/A 3.6P + 2.4 4.2RP + 2.4 16R + 2.4 
/6 4.77P + 2.4 2.7 RP + 2.4 8.5 RP + 2.4 


studying of the influences of adaptation of surface 
colours. The standard error of a single match was 
about two hue steps. The influence of small diffe- 
rences in the colour perception of the two eyes can be 
overcome by repeating the matching with the other 
eye. This repetition should also reduce the effect of 
observational errors by about 30 per cent. 

The experiment also provides evidence of the 
power of the statistical approach to experimental 
work in this field. The value of designing and analy- 
sing colour experiments according to the factorial 
scheme can be readily assessed from the experiment : 
With less than three hundred matches it was possible 
to sort out the contribution of six factors and eigh- 
teen important interactions, and to evaluate the in- 
fluences of several factors with standard errors less 
than one hue step. Where interactions were involved 
the estimation of the influences of the interacting 
factors was inevitably less accurate ; but these inter- 
actions would not have been discovered had the expe- 
riment been otherwise designed. 

The numerical results given above also serve to 
indicate the number of repetitions which would be 
required to reach a specified degree of precision in 
detailed experimentation with surface colours to help 
resolve some of the intricate problems of colour vi- 
sion. 
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Spherical Aberration of the Eye 
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Ivanorr’s data on the spherical aberration of the 
eye [1, 2] have been criticized by Koomen, Tousey, 
and ScoLnik [3]. The disagreement between IvANorr’s 
data and other determinations lies largely in what each 
experimenter considers to be the point conjugate to 
the fovea with respect to the central zone of the pupil. 
KoomEN et al. do not give any value for the central 
zone While IvaNnorr selects the zero point as represen- 
ting the refractive error of the central axial rays. 
IVANOFF’s apparatus actually does not allow him to 
measure the position of the point conjugate to the 
fovea with respect to the narrowest paraxial zone. 
The particular shapes of his curves have their basis, 
therefore, in the assumption that the accommodative 
mechanism of the eye is always so adjusted that, so 
far as the central zone of the pupil is concerned, per- 
fect conjugacy exists between the object and the 
fovea under all the conditions of accommodative sti- 
mulation used by him. 
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1953 measurements of the spherical aberration 
of the eye [2], average of the horizontal meri- 
dian of 10 eyes for accommodative stimulation 
of 0,1.50 and 3.00 diopters, The dotted lines repre- 
sent IvaNoFe’s interpretation of the data and the 
solid lines the author's interpretation. 


It is in fact impossible to determine by ordinary 
means just how the refractive state of the eye changes 
in the narrowest paraxial zone. The best estimate 
that can be made is to find the point of intersection in 
object space of two rays, one passing through the cen- 
ter of the pupil and the other coming in just a little 
eccentrically, which the subject reports or which are 
otherwise determined to be intersecting in the fovea. 
Ivanorr’s data points for 0.5 mm eccentricity actually 
constitute such measurements and should be regar- 
ded as the best approximation to the refractive condi- 
tion of the central zone of the pupil. 

When Ivanorr’s 1953 data are redrawn on the basis 
of this analysis, the following diagram is obtained 
(fig. 1). This shows that the actual change in spherical 
aberration during accommodation is small. 

It is necessary to say a word about the spacing of 
the lines on the graph. The line corresponding to zero 
accommodative stimulation indicates that on the ave- 
rage the eye is nearly 0.75D myopic. For 1.50D stimu- 
lus the accommodative response is Just about correct 
while for a 3.00D stimulus the eye remains a little 
over 0.50D hyperopic, i. e., underaccommodated. 
This finding is compatible with the observation [4] 
that most subjects will show a small amount of accom- 
modation when they are presented with a viewing 
situation at optical but not physical infinity. Simi- 
larly, most situations requiring a high accommodative 
effort are responded to with somewhat less of an ac- 
commodative effort than theoretically necessary. 

Preliminary experiments by the present writer sup- 
port the above view. The subjective optometer device 
designed by Fry [5] was used which measures in 
effect the intersection in object space of two parallel 
sheets of rays intersecting in the retina. The separation 
of the sheets in the plane of the pupil was kept cons- 
tant at§1.5 mm. Moving the pupil back and forth 
across the incident beams, no significant change in 
refractive power of the eye was found until the beams 
came near the edge of the pupil, and increasing the 
level of the accommodative response, the necessary 
adjustment in the setting of the optometer having 
been made, did not change these findings. The rela- 
tionship between accommodative stimulation and 
response, incidentally, was just as expected. There 
was a small amount of overaccommodation at low 
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levels of accommodative stimulation and an underac- 
commodation of up to 1.00D at high levels of accom- 
modative stimulation. With a separation of 1.5 mm 
between the two beams one would expect some indica- 
tion of the invagination of the pupil zone versus refrac- 
tive power curves shown by Ivanorr. None, however, 
was observed. 
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Information 


Provisional arrangements for 1956 Meeting of International Commission for Optics and accompanying 


optical symposium on “ Frontiers in Physical Optics ’’. 


Sponsors : National Academy of Sciences-National 
Research Council of the USA, National Science Foun- 
dation of the USA, American Academy of Arts and 
Sciences, Optical Society of America, International 
Union of Pure and Applied Physics, and Interna- 
tional Commission for Optics. 

Organizing Committee : The planning and conduc- 
ting of the scientific meetings is the responsibility 
of the USA National Committee of the International 
Commission for Optics. 

Further particulars can be obtained from Dr S. S. 
Batiarp, Visibility Laboratory, Scripps Institution 
of Oceanography, University of California, San Diego 
52, California, U.S. A. 

Date : Wednesday, 28th March-Tuesday, 3rd April, 
1956, 

Place : Cambridge, Massachusetts. Meetings will be 
held in the Faculty Club and the Kresge Auditorium 
and Little Theatre of the Massachusetts Institute of 
Technology. 

Programme : One day will be allocated to the mee- 
ting of the International Optical Commission, with 
opportunity for a further meeting if necessary. 


Four days will be allocated to a symposium on 
« Frontiers in Physical Optics », with the following 
selected topics as the main subjects for discussion : 

1. New phase contrast and interference microsco- 


pic and measuring devices. — 2. Thin films and fil- 
ters. — 3. Meteorological optics. — 4. Nuclear op- 
tics. 


One day will be set aside for organized visits to uni- 
versity and industrial laboratories in the Boston- 
Cambridge area. 

Speakers : Invitations will be issued by the USA 
National Committee to experts in various countries 
to present main papers on each of above-noted sub- 
jects and perhaps prepared discussion also. It is, howe- 
ver, the intention to leave ample time for free discus- 
sion in which anyone can take part. 


Language : Invited papers and prepared discussion 
will be in English. 


Note : The regular spring meeting of the Optical 
Society of America will be held on 5th-7th April in 
Philadelphia, Pennsylvania, and all interested per- 
sons are invited to attend. Information on this mee- 
ting can be obtained also from Professor BALLARD. 
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